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ABSTRACT 

We investigate the suggestion that the frequencies of the kilohertz and other QPOs in 
the emission from neutron stars in low-mass X-ray binaries are generated by geodesic 
motion of gas clumps around the star. First we assume, following previous work, that 
the dominant frequencies produced by such motions are the azimuthal frequency i^j^ , 
the apsidal precession frequency :^^p, and the first overtone of the nodal precession 
frequency v^^- We investigate whether geodesies can be found for which these fre- 
quencies agree with the observed frequencies. Correcting calculational errors made 
when the geodesic precession hypothesis was first proposed, we find that v^^ and J^^p 
disagree qualitatively with the frequencies of the kilohertz QPOs for infinitesimally or 
even moderately eccentric geodesies. These frequencies are similar to the frequencies 
of the kilohertz QPOs only for highly eccentric geodesies, with apastron to perias- 
tron ratios ^3-4; for these geodesies, 2iy-^-p differs qualitatively from the frequencies of 
the low-frequency QPOs. Next we investigate whether these frequencies would be the 
dominant frequencies produced by orbiting clumps. We find that they are not and that 
the dominant frequencies are instead harmonics and sidebands of the radial epicyclic 
frequency v,-. Finally, we show that gas dynamical constraints restrict the radii of or- 
biting clumps to be ^ 10^'^ of the stellar radius; the fractional modulation that could 
be produced by clumps is therefore -C 10^^. We conclude that there are significant 
difficulties in interpreting the kilohertz and low-frequency QPOs as consequences of 
geodesic motion of gas clumps. 
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1 INTRODUCTION 

Observations of low mass X-ray binaries (LMXBs) using the 
Rossi X-ray Timing Explorer (RXTE) satellite have led to 
the discovery of prominent, narrow features in power spec- 
tra of the accretion-powered X-ray emission of these sources, 
with frequencies ranging from ~5 Hz to ~1300 Hz (for a re- 
view, see van der Klis 2000). High-frequency quasi-periodic 
oscillations (QPOs) have been observed with frequencies 
ranging from ~500 Hz to ~1300 Hz. These 'kilohertz QPOs' 
often appear as a pair of peaks with frequencies u\ and U2 
(> vi) that move together and can shift up and down in 
frequency by up to a factor of two within a few hundred 
seconds, probably because of changes in the accretion rate. 
As they shift, v\ and U2 follow a narrow track in the V\-U2 
plane that is fixed for a given source. The kilohertz QPOs 
have frequencies in the range expected for orbital motion 
near neutron stars and are very likely a strong-field general- 
relativistic phenomenon. 
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In six kilohertz QPO sources, oscillations with frequen- 
cies of hundreds of Hertz have been observed during X-ray 
bursts. These 'burst oscillations' are almost certainly caused 
by rotation of brighter regions of the stellar surface at the 
spin frequency of the star (Strohmayer et al. 1996, 1998; 
Bildsten 1998; Miller 1999; Strohmayer 1999). Observations 
indicate that the principal burst oscillation frequency u-^^ is 
the stellar spin frequency v^ or its first overtone (see, e.g., 
Strohmayer & Markwardt 1999; Miller 1999). 

In four sources, burst oscillations and two simultaneous 
kilohertz QPOs have both been detected with high confi- 
dence. In 4U 1702-43 (Markwardt, Strohmayer & Swank 
1998), 4U 1728-34 (Strohmayer et al. 1996; Mendez & van 
der Klis 1999), and 4U 1636-53 (Zhang et al. 1996; Wij- 
nands et al. 1997a), the difference A// = !^2 — i^i between the 
frequencies of the two kilohertz QPOs is equal (to within 
4-15%) to v^^. In KS 1731-26 (Wijnands & van der Klis 
1997), Av is equal (to within <, 0.7%) to one-half u-^^. Careful 
analysis has shown that in 4U 1636—536, Av is slightly but 
significantly smaller than the spin frequency inferred from 
its burst oscillations (Mendez, van der Klis & van Paradijs 
1998). 
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In some kilohertz QPO sources, the separation fre- 
quency Au decreases systematically by 30-100 Hz, depend- 
ing on the source, as 1^2 increases by a much larger amount 
(Sco X-1: van der Klis et al. 1997; 4U 1608-52: Mendez et 
al. 1998; 4U 1735-44: Ford et al. 1998; 4U 1728-34: Mendez 
& van der Klis 1999; see also Psaltis et al. 1998). 

The fact that v^^ is close to one or two times Av in 
four sources as well as other evidence motivated the de- 
velopment of the sonic-point beat-frequency (SPBF) model 
(Miller, Lamb & Psaltis 1998; Lamb & Miller 2000). In this 
model the frequency V2 of the upper kilohertz QPO is close 
to but less than the general-relativistic orbital frequency j/^^,,^ 
at the sonic radius, where the flow in the disc changes from 
nearly circular to rapidly inspiraling, whereas the frequency 
ui of the lower kilohertz QPO is comparable to but less than 
the beat frequency i^b = t^^^^, — v^ of the neutron star's spin 
frequency with v^^^^ . The transition to hypersonic radial in- 
flow that occurs at the sonic radius Rs in this model is a 
strong- field general relativistic effect. In the SPBF model, 
the low-frequency QPOs seen in the Z and atoll sources are 
produced by interaction of the stellar magnetic field with 
the disc flow well outside the sonic point, via the magneto- 
spheric beat-frequency mechanism (Alpar & Shaham 1985; 
Lamb et al. 1985). 

The SPBF model explains naturally the existence of 
only two principal kilohertz QPOs in a given source, the 
near commensurability oi Av and i^^ , and the high frequen- 
cies, coherence, and amplitudes of these QPOs. The SPBF 
model predicts that Av is less than but close to u^ and ex- 
plains naturally the decrease in Av with increasing U2 seen 
in some sources (Lamb & Miller 2000). The inward drift of 
the accretion flow near Rb causes ui to increase faster than 
V2 as the accretion rate increases. Depending on the source, 
a differential increase of 7% is sufficient to explain the largest 
observed decrease in A//. 

A strong prediction of the SPBF model is that the neu- 
tron star in 4U 1636—536 is not spinning at the ~580 Hz 
frequency of the principal burst oscillation in this source but 
instead at ~290 Hz. This prediction ran counter to what 
was widely believed at the time (see, e.g., Stella & Vietri 
1999), but was confirmed by detection of a ~290 Hz sub- 
harmonic with an amplitude equal to 40% of the amplitude 
of the ~580 Hz oscillation (Miller 1999). The SPBF model 
also predicts that the V2-M relation may flatten in some 
sources when the region of nearly circular flow reaches the 
innermost stable circular orbit (ISCO). This signature of the 
ISCO may have subsequently been observed in 4U 1820—30 
(Zhang et al. 1998; Kaaret et al. 1999). In the SPBF model, 
the kilohertz QPOs are a strong-field general relativistic ef- 
fect and are therefor sensitive probes of the properties of the 
spacetime near the neutron star, including whether there is 
an ISCO and the gravitomagnetic torque produced by the 
spin of the star. Hence, if the model is shown to be correct, 
measurements of the kilohertz QPOs can be used not only to 
determine the properties of neutron stars but also to explore 
gravitational effects in the strong-field regime. 

Following the development of the SPBF model, Stella & 
Vietri (1998; see also Morsink & Stella 1999) proposed that 
the 10-60 Hz low-frequency QPOs (LFQPOs) and bumps 
observed in power spectra of the X-ray emission from the 
atoll sources and, possibly, the horizontal branch oscilla- 
tions (HBOs) observed in the Z sources (which are observed 



Table 1. Proposed QPO frequency identifications. 



QPO Frequency Geodesic Frequency 

'^HBO or '^L Nodal precession frequency u-^-p or "iv-^-p 

v\ Apsidal precession frequency Vj^^, 

V2 Azimuthal frequency v-^ 



simultaneously with the kilohertz QPOs; see van der Klis 
2000) might be the first or second harmonics of the nodal 
precession frequency v^^, of gas clumps moving on slightly 
inclined, infinitesimally eccentric geodesies (lEGs) with az- 
imuthal frequencies equal to 1/2, around neutron stars with 
the 250-350 Hz spin rates inferred from the burst oscilla- 
tions and the SPBF model. For such stars, nodal precession 
is driven by the prograde gravitomagnetic (Lense-Thirring) 
torque, but is partially offset by the retrograde torque pro- 
duced by the star's rot at ion- induced quadrupolar distortion. 

Recent work on the global modes of viscous accretion 
discs (Markovic & Lamb 1998) has demonstrated that such 
discs have a dense spectrum of gravitomagnetically precess- 
ing modes that are localised near the inner boundary of the 
region of nearly circular flow. The highest-frequency of these 
modes precesses with a frequency only slightly lower than 
the nodal precession frequency u^p of a test particle orbit- 
ing at the same radius. Contrary to what had been expected 
since the pioneering work of Bardeen & Petterson (1975), the 
first few of these modes are very weakly damped and may 
therefore be excited, producing detectable variations in the 
X-ray emission (Markovic & Lamb 1998; see also Armitage 
& Natarajan 1999). However, further work (Psaltis et al. 
1999a; Morsink & Stella 1999; Kalogera & Psaltis 1999) has 
shown that the frequencies and frequency behaviors of the 
LFQPOs and HBOs are inconsistent with the proposal of 
Stella & Vietri (1998) that these frequencies are the first or 
second harmonics of the nodal precession frequency of the 
lEG with azimuthal frequencies ~ V2 around stars with the 
250-350 Hz spin rates inferred from the burst oscillations 
and the SPBF model. 

Stella & Vietri (1999; hereafter SV) subsequently pro- 
posed that i/2 and fi are the azimuthal and apsidal pre- 
cession frequencies v^^ and fj^p of gas clumps orbiting the 



star on eccentric geodesies; Vj^ 



fr in terms of the 



radial epicyclic frequency fr. They again proposed that the 
frequencies i^hbo and v^ of the low-frequency QPOs and 
bumps seen in power spectra of the Z and atoll sources are 
the first or second harmonics of the nodal precession fre- 
quency Vj^p of these same geodesies, but argued that u^ is 
may be S> 300 Hz to give better agreement of u^p with !/hbo 
and i^L- The frequency identifications proposed by SV are 
listed in Table Ol. SV simply assumed that radiating or ob- 
scuring gas clumps with the required properties form near 
the star; that hydrodynamic, radiation, and magnetic forces 
are negligible, so that the gas clumps move on geodesies; 
that 2iy^p , Vj^p , and i/^ are the dominant frequencies gener- 
ated by clumps in purely geodesic motion about a neutron 
star; and that all the clumps move on geodesies with the 
same periastron and apastron radii rp and r^- The large 
variations of i^hbo or z^^i '^ij and V2 that occur in the indi- 
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vidual sources are attributed to gas clumps being fed onto 
different geodesies at different times. 

SV first pointed out that the frequencies z/^p and i/^^ of 
infinitesimally eccentric geodesies (lEGs) around nonrotat- 
ing stars with masses of 1.8-2.2 M0 are comparable to the 
values of i^i and 1^2 observed in Sco X-1, although they do 
not track the tightly correlated variation of vi and U2 (see 
their Fig. 1). Stella & Vietri then argued that sequences of 
moderately eccentric geodesies (MEGs) would give apsidal 
precession and azimuthal frequencies more consistent with 
the observed frequencies of the kilohertz QPOs. 

Once geodesies of finite eccentricity are allowed, the 
geodesic precession hypothesis does not by itself predict a 
relation between Uj^^ and u^ . The reason is that the charac- 
teristic frequencies of such geodesies are functions of two 
parameters: rp and ra. The set of geodesies that satisfy 
ra. > rp > Re generate frequencies that cover a wide area 
in the two-dimensional i^j^p-v^ plane, rather than a one- 
dimensional curve. In order to obtain a z^^p'^^k relation that 
could explain the tightly correlated variation of i/i and V2, 
some physical constraint is needed that will convert the two- 
parameter set of allowed geodesies to a one-parameter al- 
lowed sequence of geodesies that generates a curve in the 
^Ap"^K plane relating z/^p and v^^. 

SV addressed this issue by proposing that some physi- 
cal mechanism keeps rp constant in a given source while ra 
varies. They suggested that rp may remain constant because 
it is equal to the radius risco of the innermost stable circular 
orbit, to the radius Vm at which the magnetic field of the star 
first couples to the accreting gas, or to the equatorial radius 
Rb of the neutron star. However, the periastron radius of an 
eccentric geodesic can be significantly smaller than risco (see 
Markovic 2000) and the magnetic coupling radius Vm varies 
as the accretion rate varies (see Lamb 1989). 

SV reported that they were able to construct MEG se- 
quences with apsidal precession and azimuthal frequencies 
that agree approximately with the kilohertz QPO frequen- 
cies observed in Sco X-1, assuming that rp is constant but 
can have a ny v alue greater than risco- However, as we explain 



in Section 3.1 and the Appendix, SV computed ly^^ incor- 
rectly. Karas (1999) noticed this error (see his footnote 2), 
but thought that its effect was small because the geodesies 
that SV considered have small eccentricities e ~ 0.1, which 
corresponds to ra/rp — 1.2. However, in order to fit the data, 
Vr must be ;$0.3 2/2 • A given clump therefore makes only a 
small fraction of its full radial excursion during each orbit 
of the star and hence — in contrast to the behavior of parti- 
cles in similar Newtonian orbits — the variation with radius 
of its rate of azimuthal phase advance does not average out. 
The error made in SV therefore produced a 10% (« 100 Hz) 
error in v^^ , which is very large compared to the uncertain- 
ties in 1/2 for Sco X-1. When the correct values of i/^ are 
used, the MEG sequences cited by SV give Vj^p-i/^ relations 
that disagree qualitatively with the Sco X-1 frequency data 
(see Section 3.1). SV did not consider the frequency data 
available on other sources. 

Motivated by the proposal of SV, Karas (1999) investi- 
gated whether geodesic sequences can be constructed which 
give i-'r-i^K relations that agree with the A:y-i/2 correlations 
observed in Sco X-1 and 4U 1608—52, if M is treated as 
a free parameter for each source and rp and ra are freely 
chosen for each measurement of vi and i-'2. Karas also in- 



vestigated whether geodesic sequences can be constructed 
which give v^-Vj^^ relations that agree with the l\v-U2 corre- 
lations observed. This assumes that the principal kilohertz 
QPO frequencies are Uy^ — 2ur and v^ ~ 1^1', its physical 
motivation is unclear. Karas considered only nonrotating 
masses (he did not construct any stellar models) and allowed 
highly eccentric geodesies (e up to 0.4, which corresponds to 
ra./rp — 2.3). Not surprisingly, given that varying rp and 
ra allows coverage of a large area of the i^r-J^K plane (see 
above), he found that fits are possible, although he did not 
report the parameters or the value of x^/dof for any indi- 
vidual fits. Such fits imply that ra is always a function of 
rp, i.e., that specifying rp uniquely specifies a particular ec- 
centric geodesic. Moreover, the function ra = ra(rp) must 
be different for every source. The physical motivation for 
this is unclear. For some fits, Karas assumed that a clump 
only needs to complete 92%-97% of the radial extent of its 
geodesic, because of the clump's finite size. The physical 
basis of this assumption is unclear. 

Psaltis, Belloni & van der Klis (1999b) have pointed to a 
possible correspondence between features in the power spec- 
tra of black holes and neutron stars in LMXBs, emphasising 
the similar properties of the 0.1-100 Hz QPOs observed in 
the atoll sources, Cir X-1, and some black hole sources and 
the 15-60 Hz HBOs seen in the Z sources. Somewhat more 
speculatively, Psaltis et al. suggested that the power spec- 
tral 'humps' sometimes observed at 10-60 Hz in black hole 
sources may be related to the lower-frequency members of 
the two kilohertz QPOs seen in the Z and atoll sources, 
which have frequencies vi ~ 300-800 Hz. 

Motivated by the results of Psaltis et al. (1999b), Stella, 
Vietri & Morsink (1999; hereafter SVM) compared the fre- 
quencies of lEGs around neutron stars and black holes with 
the frequencies of QPOs and other features observed in neu- 
tron star and black hole LMXBs. Assuming that u^, Vj^p, 
and 2i/^p are the dominant frequencies produced by gas 
clumps moving on such lEGs, that the spin rates of the neu- 
tron stars in LMXBs range from 300 to 900 Hz, and that the 
dimensionless angular momenta j = cJ/GM^ of the black 
holes in LMXBs range from 0.1 to 0.3 (here J and M are 
the angular momentum and mass of the black hole), SVM 
showed that there are lEGs around 1.95 M© neutron stars 
which give values of v^, i^^pi and 2Vj^p that agree within 
factors ~2 with the frequencies of the kilohertz and lower- 
frequency QPOs observed in neutron star LMXBs and that 
there are lEGs around black holes which give values of 2u-f^p 
and Uj^p that agree within factors ~2 with the frequencies of 
the 0.1-10 Hz and 2-300 Hz features observed in the power 
spectra of some black hole LMXBs. However, SVM did not 
attempt to fit any frequency data. 

Recently, Psaltis & Norman (2000) proposed a model 
of QPOs in neutron star and black hole LMXBs based on 
the assumption that there is a "resonant ring" in the inner 
accretion disc. A ring with the properties proposed has a 
dense spectrum of resonances, including resonances at the 
frequencies 2v^ 



NP ' "^AP ' 



and v-y of an lEG at the radius of the 
ring. Psaltis & Norman assumed that these three frequen- 
cies are the dominant frequencies produced by the resonant 
ring and identified them with the frequencies of the low- 
frequency and kilohertz QPOs in the neutron star LMXBs. 
If it could be established that the frequencies of the 
most prominent QPOs and other features observed in power 
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Table 2. Types of the geodesies considered in the present work. 



Infinitesimally Eccentric Geodesies (lEGs) Geodesies with Finite Eccentricities" 



Quantity 



Periastron radius rp Equal to r^ 



Apastron radius ra 



Freely chosen for each frequency pair 



Fixed for a given source; varies from source to source 
Freely chosen for each frequency pair 



"We refer to geodesies with finite eccentricities as moderately eccentric geodesies (MEGs) if ra < 1.5 r'p (e < 0.2) or as highly 
eccentric geodesies (HEGs) if ra > 1.5 rp (e > 0.2). 



spectra of the X-ray emission from neutron star and black 
hole LMXBs are simply the characteristic frequencies of test 
particle geodesies around these objects, measurements of 
these frequencies would provide a very simple way to ex- 
plore the effects of strong-field gravity. Moreover, as dis- 
cussed by Markovic (2000), the periastron radius of an ec- 
centric geodesic can be smaller than the radius of the ISCO. 
Hence, if the highest frequency QPOs could be shown to be 
produced by clumps of gas moving on geodesies with sub- 
stantial eccentricities, it might be possible to derive even 
tighter upper bounds on the radii of the neutron stars in 
the kilohertz QPO sources and hence on the hardness of the 
equation of state of neutron star matter than have been de- 
rived from the models in which the frequencies of the upper 
kilohertz QPOs are the orbital frequencies of gas in nearly 
circular orbits (Miller et al. 1998; Schaab & Weigel 1999). 
On the other hand, by rejecting any causal connection be- 
tween the frequency difference Ai/ and the stellar spin fre- 
quency v^ , such a model would leave unexplained the origin 
of the burst oscillations and the closeness of Av to the burst 
oscillation frequency. 

In the present paper we investigate further the possibil- 
ity that the frequencies of various QPOs and other features 
in power spectra of the X-ray emission from neutron stars in 
LMXBs are certain frequencies of special geodesies around 
these stars. For the purposes of this analysis, we simply as- 
sume that gas clumps with the required properties can be 
formed and destroyed in the required way; that they are in- 
jected onto the required geodesies; and that once injected, 
the motion of the clumps is unaffected by radiation, mag- 
netic, or gas pressure forces, so that their motion is purely 
geodesic. We then explore whether geodesic frequencies are 
consistent with the observed QPO frequencies; compute the 
oscillation frequencies that would be produced by emitting, 
reflecting, or obscuring clumps moving on geodesies around 
a neutron star; and analyze the constraints imposed on such 
models by basic gas dynamics. We consider the geodesies 
and frequency identifications proposed by SV, SVM, Karas 
(1999), and Psaltis & Norman (2000), but we do not restrict 
ourselves to these models. The various types of geodesies we 
study are listed in Table g. 

In Section Q we describe how we compute the neutron 
star models, spacetimes, and geodesies that we need and 
the procedure we use to construct the geodesic frequency 
relations that fit best the measured QPO frequencies. 

In Section a we investigate whether certain geodesic fre- 
quencies are consistent with the frequencies of the power 
spectral features observed in the atoll and Z sources. For 
the purposes of this analysis, we simply assume that clumps 
moving on geodesies produce X-ray oscillations with the fre- 
quencies f^ , Uf^-p , and 2i/j^p , and then examine whether one 



can construct geodesies for which these frequencies agree 



with the observed frequencies. We find that the f^ 



rela- 



tions given by the best-fitting sequences of lEGs and MEGs 
are qualitatively different from the observed iyi-1/2 corre- 
lations, for any acceptable neutron star models and spin 
rates. We therefore explore whether dropping the require- 
ment that geodesies be only moderately eccentric would al- 



low acceptable fits. We show that i^^ 



relations similar 



to the observed kilohertz QPO frequency correlations are 
possible using highly eccentric geodesies (HEGs), but only 
if the geodesic sequences include HEGs with apastron to 
periastron ratios > 3; most of the best-fitting frequency re- 
lations are, however, formally unacceptable. We find that 
frequencies other than i^^p and i^^ produce frequency re- 
lations that disagree qualitatively with observed kilohertz 
QPO frequency correlations. 

Next, we examine whether sequences of HEGs can be 



found that give 2i/j., 



and ly. 



relations that are si- 



multaneously consistent with the frequencies of the low- 
frequency QPOs and the kilohertz QPOs. We find that the 



best-fitting 2;/^^ 



and i^AP'^^K relations disagree qualita- 



tively with the observed correlations, for realistic neutron 
star models and physically allowed spin frequencies. 

We also investigate whether there are sequences of 



geodesies that give 2Uj^ 



relations consistent with the 



frequencies of the low- and high-frequency spectral features 
observed in Cir X-1. Geodesies can be found that give rela- 
tions consistent with most of the Cir X-1 frequency data, but 
only geodesies around neutron stars with very low masses 
(M ~ 1.3 Mq) and high spin rates (900 Hz) have frequencies 
that also approach the frequencies of the highest-frequency 
QPOs observed in Cir X-1. 

In Section H we compute numerically the X-ray oscilla- 
tion frequencies that would be generated by occulting, re- 
flecting, and radiating clumps moving on geodesies around a 
neutron star by simulating the waveforms produced by such 
clumps. We flnd that contrary to what has been assumed 
in previous work, most of the power generated by orbiting 
clumps appears in peaks of roughly equal total power at the 
radial epicyclic frequency fr, at i/^ -\- fr, and at v-^ -\- 2i/r. 
There is power at v^ , but it is always <^ 25% of the power 
in each of the three dominant peaks and there are often 
many other peaks as strong or stronger than the peak at 
i/j^ . For the relevant geodesies, the power at Vj^^ and at v^p 
generated by most mechanisms is <, 10% of the power in 
the dominant peaks. We were unable to find any effect or 
geometry that produces significant power at 2!/„p . 

In Section |5| we explore the constraints on orbiting 
clump models and on the amplitude of the X-ray modulation 
that can be generated in such models imposed by basic gas 
dynamics. We find that gas dynamical constraints require 
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the size of an individual clump to be a fraction < 10"'' of 
the stellar radius and the density of the gas in the clumps 
to be 3> 10 times the interclump density and S> 10 times 
the mean density near the star. As a result, the maximum 
X-ray modulation that clumps could produce is ^ 10"". 

In Section M we discuss our principal results and con- 
clude that there are significant difficulties with the geodesic 
precession and resonant ring models of the kilohertz and 
other QPOs. 



2 COMPUTATION OF GEODESICS AND 

FREQUENCY COMPARISON PROCEDURE 

2.1 Computation of geodesies 

We first computed numerically a large suite of neutron star 
models and spacetimes for a range of neutron-star matter 
equations of state (EOS), stellar spin rates, and masses. We 
then integrated the equations of bound geodesic motion (see 
Markovic 2000) in the exterior spacetimes of these neutron 
star models to determine the characteristic frequencies and 
other properties of possibly relevant geodesies. 

Modern numerical codes (Komatsu, Eriguchi & Hachisu 
1989a, 1989b; Cook, Shapiro & Teukolsky 1992, 1994a, 
1994b; Stergioulas & Friedman 1995; Nozawa et al. 1998) al- 
low accurate construction of uniformly spinning, equilibrium 
stellar models using any tabulated EOS. We use the partic- 
ular variant developed by Stergioulas & Friedman (1995). 
With the exception of Cir X-1, the neutron star models that 
give the best fits to the frequency correlations we consider in 
this paper all have masses > 1.7 Mq and spin rates < 600 Hz, 
substantially lower than their mass-shedding limits. 

We considered a range of neutron star matter equa- 
tions of state. We find that stars constructed using the mod- 
ern, realistic EOS A18-|-UIX-|-5iib (Akmal, Pandharipande & 
Ravenhall 1998; Pandharipande, Akmal & Ravenhall 1998) 
give geodesies that are almost identical to those given by 
the earlier realistic UU EOS of Wiringa, Fiks & Fabrocini 
(1988). In order to facilitate comparison with previous work, 
we mostly report results based on the UU EOS. In order 
to explore the effects of uncertainties in the EOS, we also 
studied stars constructed using the early EOS of Bethe & 
Johnson (1974; EOS C in the compilation by Arnett & Bow- 
ers 1977), which is softer than the UU EOS at high densi- 
ties. In order to allow comparison with previous work that 
used this EOS (e.g., SV), we also investigated stars based on 
an early mean-field EOS constructed by Pandharipande & 
Smith (1975; EOS L in Arnett & Bowers 1977), even though 
this EOS is now considered unrealistically hard. These three 
equations of state are discussed in more detail in the com- 
panion paper (Markovic 2000) . 



2.2 Construction of frequency relations 

The frequencies i^np) ^api ^-^d ^'k ^^ ^ geodesic are func- 
tions of its periastron and apastron radii rp and ra (we use 
circumferential radii everywhere), the EOS of neutron star 
matter, and the mass M and dimensionless angular momen- 
tum j = cJ/GM^ of the neutron star. For a given EOS and 
M, j is a function of the neutron star spin rate Vs, which 



may be directly observable. We therefore parameterise the 
effects of stellar rotation by Us . 

For infinitesimally eccentric geodesies, we assume that 
the sequence {ra,i} of apastron radii is the same as the se- 
quence {^p^i} of periastron radii and then use the maxi- 
mum likelihood method to determine the {rp,i} sequence 
that gives the z^ap'^k sequence that best fits the observed 
i/i-z^2 sequence. We treat each rp,i as a free parameter, re- 
quiring only that it exceed risco. 

For geodesies with finite eccentricities, the ith mem- 
bers of the sequences {i^Np;}, {i^apz}' ^^'^ {^k t} depend 
on both rp,i and ra,i. Hence the geodesic precession hypoth- 
esis does not by itself predict a relation between the different 
characteristic frequencies of a geodesic. As explained in the 
Introduction, some physical constraint is needed that will 
convert the two-parameter set of geodesies specified by rp 
and Ta to a single-parameter sequence of geodesies that will 
generate a curve in, e.g., the z^af'^^k plane. We follow Stella 
& Vietri (1999) in assuming that rp is constant in a given 
source but that it can be chosen freely for each source to 



give the best possible fit of the fci^j., 



and v. 



rela- 



tions to the observed frequency correlations (we require only 
that rp exceed the equatorial radius Re of the star, if there 
is no marginally bound geodesic, or that rp exceed r^h, the 
periastron radius of the marginally bound geodesic). Here 
A; (= 1 or 2) is the harmonic of the nodal precession fre- 
quency that the hypothesis identifies with ^'hbo or Uj^. For 
a given rp, the frequencies i^^pj ^api ^^^ ^k ^^^ continuous 
functions of ra. We assume that the ra sequence for each 
source can be chosen freely to give the best possible agree- 
ment of the kVj^p, j/^p, and u^ sequences with, e.g., the 
observed i^bo, t^i, and V2 sequences. We determine which 
ra sequence agrees best with the frequency data using the 
maximum likelihood method. We refer to geodesies with fi- 
nite eccentricities as moderately eccentric if ra,max/rp < 1.5 
and as highly eccentric if ra,max/rp > 1.5 (see Table bl). 

For each EOS, we treat M as a free parameter within 
the mass range allowed by the EOS. The variation of i^^p 
with i/j^ is relatively insensitive to the spin rate of the star, 
unless the stellar mass required is near the maximum stable 
mass, in which case nonzero spin makes the fits significantly 
worse. In contrast, the variation of v^p with v-^ is sensitive to 
the spin of the star. In order to give the geodesic precession 
hypothesis the best possible chance of fitting the data, we 
treat Vs as a free parameter, allowing any value between 
zero and the mass-shedding limit, even though the neutron 
stars in the kilohertz QPO sources are thought to have spin 
frequencies ~250-350 Hz. 

We assume that a necessary (although clearly not suffi- 
cient) condition for the geodesic precession hypothesis to be 
acceptable is that the best-fitting sequence of geodesies gives 
frequencies consistent with the QPO frequencies observed in 
all the sources to which the hypothesis is supposed to apply. 
In assessing whether the frequencies given by a particular 
geodesic sequence are consistent with the observed frequen- 
cies of the relevant QPOs, we use the frequency uncertainties 
reported in the literature. These uncertainties generally in- 
clude only statistical errors. We caution that when a QPO 
or other feature is weak, systematic errors introduced, for 
example, by uncertainties in the continuum, are likely to be 
significant. 
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2.3 Fits to kilohertz QPOs 

Consider first the kilohertz QPOs. Suppose there are n si- 
multaneous measurements of the two kilohertz QPO fre- 
quencies vi and 1^2 and let vi^i and !/2,i denote, respectively, 
the ith (simultaneous) pair of measurements. We assume 
that the errors in Vi^i and i/2,i are normally distributed with 
standard deviations ai^i and a2,i- Then the likelihood of any 
given set of the parameters M, v^, {rp,i}, and {ra,i} is 



/;(M,i/^,{rp,,},{ra,,}) = - 



n 



2-11 (J li (J-ZA 



X exp 



[viA-V^F^if {l'2,: 



M,. 



^<^h 



(1) 



where {Vj^^ .} and {i'k ;} are the apsidal precession and az- 
imuthal frequency sequences given by the EOS, M, v^, and 
the radius sequences {rp^i}, and {ra,i}. 

We characterise the goodness of a fit by x^/dof = ji^ /m, 
where x^ is twice the sum of the squared terms in the ex- 
ponential in equation (hi) and m is the number of degrees of 
freedom. For n data points (2n frequencies), fits of infinites- 
imally eccentric geodesies have n-\-2 parameters (the {ra,i} 
plus M and u^), or n + 1 parameters if the spin is fixed a 
priori, so m = 2n — n — 2 = n — 2, or n — 1 if the spin is 
fixed; fits of geodesies with finite eccentricities have n -I- 3 
parameters (the {ra.i} plus rp, M, and u^) or n -1- 2 if the 
spin is fixed, so m = n — 3, or n — 2 if the spin is fixed. 

We begin by picking an EOS and spin frequency Vs . The 
frequencies u^-p and u-^ depend only on the exterior space- 
time close to the rotation equator and are single- valued, 
continuous functions of rp and r^,. As shown in the com- 
panion paper (Markovic 2000), the exterior spacetime close 
to the rotation equator of a star with mass M and spin 
frequency Vs can be described by the eight metric func- 
tions po(M, i/^,s), ao{M,u^,s), f3o{M,u^,s), 'yo{M,u^,s), 
p{M,h'^,s), a{M,u^,s), P{M,v^,s) and ^{M,v^,s), where 
s = r/{r + Ri,). We determine these metric functions for each 
v^ and each Ad < Mmai^ii's) that we consider, by interpolat- 
ing in a sequence of metric functions computed previously 
for this value of u^ and a sufficiently dense sequence {Mi} of 
stellar masses. Here Mma^ii'^) is the maximum stable mass 
for stars of spin frequency v^ . 

Once the metric functions are known, we can compute 
i/^p and i/j^ as functions of rp and ra. For infinitesimally 
eccentric geodesies, we determine numerically the sequence 
{ra,i} of Ta values that, with rp,i = ra,i, maximises jC. 
For geodesies with finite eccentricities, we determine the se- 
quence {ra,i} that maximises jC for a given rp. We then vary 
rp and M to find the values that give the largest maximum 
value of C. As noted above, the i^^p'^^K relation given by a 
particular sequence of geodesies is fairly insensitive to v^ and 
hence the quality of the fit to the 1/1-1/2 correlation observed 
in a given source varies only slightly over a wide range of 
spin rates. We therefore do not attempt to determine the 
best-fit spin frequency accurately, but instead quote a range 
of spin frequencies that have similar likelihoods. 

Although we fit the i/j^p-Vj^ relations allowed by a given 
geodesic precession hypothesis to the 1/1-1/2 data as de- 
scribed above, in our figures we compare the relation be- 
tween 1/,- — 1/^ — i/j^p and 1/^ given by the best-fitting se- 



quence of geodesies with the observed correlation between 
Ai/ = U2 — V1 and V2 , because one can assess the goodness of 
the fits much more easily by examining such plots. The un- 
certainties in /S.U shown in these plots were not used in the 
fitting procedure and are shown only to give an approximate 
visual impression of the goodness of fit; they were computed 
by adding the uncertainties in vi and U2 in quadrature. 



2.4 Fits to kilohertz and low- frequency QPOs 

Consider now the spectral features at uhbo or v-^ . We wish 
to determine whether a sequence of geodesies can be con- 



structed that gives kv-^ 



and V. 



relations that are 



simultaneously consistent with, e.g., the observed i/nBo-V2 
and U1-V2 correlations. For a given EOS, the parameters are 
the same as for the kilohertz QPOs alone, but their likeli- 
hood is now 



/:(M,!/,;{rp,,},{ra,,}) 
{Pi.i - i/^p. 



X exp 



2aiA 



J-i 2Tvaii 


0"2, 


''"HBO.i 


il/2,i - l/^,if 






2'^h 




('^HBO.i ~ ^'^NP, 


.r^ 




2'^Sbo, 






' 



(2) 



where i/^p^ = v^p{M,v^,rp^i,rs^^i). In contrast to !/^p and 
1/-^, u^p depends sensitively on u^ and hence we can deter- 
mine fairly accurately the spin frequency that gives the best 
fit. 

We again characterise the goodness of a fit by x^/dof = 
X^ /m, where x^ is twice the sum of the squared terms in 
the exponential in equation (0) and m is the number of 
degrees of freedom. For n data points (3n frequencies) , fits of 
infinitesimally eccentric geodesies have m — 2n — 2, or 2n— 1 
if the spin is fixed; fits of geodesies with finite eccentricities 
have m = 2n — 3, or 2n — 2 if the spin is fixed. 



3 FREQUENCY RELATIONS 

We now investigate whether sequences of geodesies can be 
found that have characteristic frequencies which agree with 
the frequencies of the kilohertz QPOs observed in the Z and 
atoll sources. We first explore whether it is possible to con- 
struct sequences of lEGs or MFCs (see Table bl) that give 
i/j^p-v^ relations consistent with the observed correlations of 
i/i with V2 . We find that no such sequences exist for accept- 
able neutron star models and spin rates. We also show that 
the behavior of other characteristic frequencies, such as 1/-^ 
and Vr + 1/^, is inconsistent with the observed 1/1-1/2 correla- 
tions for geodesies of any eccentricity. We therefore consider 
whether allowing highly eccentric geodesies (HEGs) would 
make it possible to construct geodesic sequences that give 
^Ap"^K relations consistent with the observed 1^1-1^2 correla- 
tions. We also explore whether there are sequences of HEGs 
with t^Ap-t'K a-nd 2i/^p-i/^ relations that are simultaneously 
consistent with the observed 1/1-1/2 and Vi^-v^ or vhbo-i/^. 
correlations. Finally, we investigate the frequency correla- 
tions observed in Cir X-1. 
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Figure 1. Upper panel: Vr-v^ relations given by two sequences 
of infinitesimally eccentric geodesies that best fit tlie v\ and V2 
data for the Z source Sco X-1. The dashed curve shows the v^-u^ 
relation given by the best-fitting geodesic sequence around a non- 
rotating star (M = 1.99 Mq; x^/dof = 38). The solid curve shows 
the relation given by the best-fitting geodesic sequence around 
a rotating star constructed using the UU EOS(M = 2.21 Mq, 
Va = 450 Hz; x^/dof = 35); the relation given by the best-fitting 
geodesic sequence around a rotating star constructed using EOS L 
is almost identical. The dotted curve shows the v^-u^ relation 
given by the best-fitting sequence of geodesies around a black hole 
(M = 5.8 M0, j = 0.89; x^/dof = 1.6). Lower panel: Circumfer- 
ential radii of the geodesies that give the frequencies plotted in 
the upper panel. 
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Figure 2. Upper panel: Comparison of the Vr-v^^ relations dis- 
cussed in the text with measurements of h.v and V2 in Sco X-1. 
The solid curve shows the t'r-i'K relation given by the geodesic 
sequence around a nonrotating star that SV reported as fitting 
the Sco X-1 data; x'^/dof for this sequence is 150. The dashed 
and dash-dotted curves show the i^r-^^K relations given by the 
geodesic sequences around the stars with spin rates of 300 Hz and 
600 Hz that SV also reported as fitting the Sco X-1 data. The 
dotted curve shows the t'r-t'K relation given by the sequence of 



geodesies around a rotating UU star with rp > 



but any 



value of ra (> rp) that best fits the vi and V2 data for Sco X-1 



(Af = 2.21 iW0, 1/^ 



450 Hz, rp 



5.52; X /dof = 21). 



Lower panel: Circumferential apastron radii of the geodesies that 
give the frequencies plotted in the upper panel. 



3.1 Frequencies of infinitesimally and moderately 
eccentric geodesies 

The most numerous and precise kilohertz QPO frequency 
measurements are those for the Z source Sco X-1 (van der 
Klis et al. 1997) . We therefore consider consistency with the 
frequency correlation observed in this source as a necessary 
condition for any QPO model to be viable (a very similar 
correlation has been observed in 4U 1608—52; see Mendez et 
al. 1998) . The frequency relations given by geodesies around 
nonrotating stars provide useful guidance because the Vi-v^ 
relation is relatively insensitive to the stellar spin rate and 
the relation for a nonrotating star doesn't depend on the 
EOS of neutron star matter, provided that it is hard enough 
to support the stellar mass being considered and gives stars 
with Re < ra.min). Only geodesies near risco are relevant, 
because Vr decreases appreciably with increasing v^^ only 
near Hsco. 

The dashed curve in Fig. |l| shows the frequency relation 
constructed from lEGs around a nonrotating star that best 
fits the Sco X-1 data; the stellar mass derived from this fit is 
1.99 Mq, which is a reasonable mass for a neutron star in an 



LMXB. The dashed curve is very similar to the curve shown 
in Fig. 1 of SV for lEGs around a nonrotating star with a 
mass of 2.0 Mq. As Fig. hi shows, the frequency relation for 
the best-fitting sequence of lEGs around a nonrotating star 
is much steeper than the observed frequency correlation; the 
qualitative disagreement between the best-fitting frequency 
relation and the frequency data is indicated by the large 
X^/dof for this fit, which is 38. 

As noted by SV, stellar rotation has little effect on the 
relation between u^^ and u^^ . This is illustrated in Fig. nl by 
the similarity of the frequency relation for the rotating star 
constructed using the UU EOS that best fits the Sco X-1 
data to the relation for the nonrotating star that fits best. 
The frequency relation for the rotating star constructed us- 
ing EOS L that best fits the data is also very similar to 
the relation for the nonrotating star that fits best. Although 
prograde spin does decrease slightly the slope of the pre- 
dicted Vr-v^ relation for a given mass, the effect is small. 
The reason is that the dragging of inertial frames causes risco 
to decrease, which significantly increases the azimuthal fre- 
quency of geodesies near Hsco for a given stellar mass. Hence 
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the mass of the star must be increased in order to keep the 
azimuthal frequencies in the observed range (all frequencies 
scale approximately as 1/M). This drives the stellar mass 
close to the maximum mass allowed by the EOS. Indeed, the 
mass of a rotating UU star that fits best is 2.21 Mq, which 
is the maximum stable mass for this EOS and the best- 
fitting spin rate (450 Hz). The x^/dof for this fit is 35.5, 
only slightly smaller than the x'^/dof for the nonrotating 
star that fits best. The fit is not any better for harder equa- 
tions of state. Although they give larger maximum masses, 
i^i falls more steeply because of their larger size and hence 
larger quadrupolar deformation (Markovic 2000). 

SV argued that i/^p and v^ would agree with the ob- 
served frequencies of the kilohertz QPOs if clumps form on 
MEGs rather than lEGs. To support their argument, they 
showed in their Fig. 2 a frequency relation supposedly for a 
sequence of MEGs with Vp = 6.25Af around a Af = 1.9 M© 
nonrotating star; this relation appears to agree fairly well 
with the Sco X-1 frequency data. However, as explained in 
detail in the Appendix, the fr-i^K relation plotted by SV was 
computed incorrectly. The solid curve in Fig. ti shows the 
actual t'r-i^K relation for the sequence of MEGs specified by 
SV. It is similar to the relation for the best-fitting sequence 
of lEGs and is much steeper than the data (x^/dof = 150). 

SV also reported that sequences of MEGs with rp ~ 
6.18M and 6.17M around rotating neutron stars with 
masses of 1.94M0 and l.QSM© constructed using EOS L 
with, respectively, u^ — 300 Hz and 600 Hz also give I'l-i'i^ re- 
lations similar to the A!/-i/2 correlation observed in Sco X-1. 
We have computed the i^r-J^K relations for these geodesic se- 
quences. The results are shown as the dashed and dashed- 
dotted curves in Fig. bt these curves lie 100-200 Hz above 
the Sco X-1 frequency correlation. 

The dotted curve in Fig. shows the frequency rela- 
tion for the best-fitting sequence of geodesies around a ro- 
tating star constructed using geodesies that meet the re- 
quirement Tp > Tisco proposed by SV but are otherwise un- 
restricted. With this freedom, the low-frequency geodesies 
become more eccentric (ra,max/''p ~ 1.5), but the frequency 
relation remains much steeper than the data (x^/dof = 21). 
For the reasons discussed above, the stellar mass that fits 
best is equal to the maximum mass for the best-fit spin 
rate. 

The frequency relations for sequences of MEGs around 
stars constructed using other equations of state are also in- 
consistent with the Sco X-1 data. The frequency relations 
given by sequences of lEGs and MEGs also disagree with 
the kilohertz QPO frequencies of other sources, although 
frequency measurements for sources other than Sco X-1 are 
sparser and less precise. For example, the best-fitting se- 
quence of MEGs with Tp > risco around a UU star with 
the 363 Hz spin rate inferred from the burst oscillations of 
4U 1728-34 (Strohmayer et al. 1996; Mendez & van der Klis 
1999) has M = 1.92 M© and gives x^/dof « 15. If instead 
the stellar spin rate is treated as a free parameter, the best- 
fitting MEG sequence with rp > risco has u^ « 900 Hz and 
M = 2.20 Mq and gives X^/dof = 9.4. This mass is again 
very close to the maximum stable mass for this EOS. 

Almost all the best-fitting neutron star masses given 
by MEG sequences are very close to the maximum stable 
mass for the assumed EOS, which is implausible. More seri- 
ously, the i^r-i^K relations given by sequences of MEGs devi- 



ate by large fractional amounts {Aur/vr as much as 30% or 
Az/j^/z/j^ as much as 15%) from the A!/-z^2 correlations ob- 
served. Any forces strong enough to bring these frequency 
relations into agreement with the data would invalidate the 
geodesic motion hypothesis. 

Fig. n] shows that one can construct a sequence of lEGs 
around a 5.8 Mq, rapidly spinning {j — 0.89) black hole that 
gives a frequency relation which agrees fairly well with the 
kilohertz QPO frequencies observed in Sco X-1 (x^/dof = 
1.6). There is compelling evidence that Sco X-1 is not a black 
hole (see, e.g., van der Klis 2000). The fairly good agreement 
of lEGs in the Kerr spacetime with the Sco X-1 frequency 
data is a reminder that successfully fitting QPO frequencies 
does not by itself validate a model. The wealth of other 
information that is available about the kilohertz QPOs and 
the kilohertz QPO sources must also be taken into acount. 

These results show that regardless of the neutron star 
model and spin rate assumed, the MEG hypothesis gives fre- 
quency relations that are qualitatively different from the ob- 
served correlations between the kilohertz QPO frequencies. 
As explaining this correlation was the primary motivation 
for introducing the MEG hypothesis (SV), we consider this 
failure serious enough to set aside the MEG hypothesis and 
we therefore do not consider it further. 



3.2 Frequencies of highly eccentric geodesies 

We consider now whether allowing highly eccentric geodesies 
(HEGs) would make it possible to construct sequences that 
give values of z/^p and v^^ consistent with the values of ui 
and 1/2 observed in the kilohertz QPO sources. As before, we 
assume that the periastron radius rp has a constant value 
> Re, if there is no marginally bound geodesic, or > rmb, 
but we allow the apastron radius ra to have any value > rp. 



3.2.1 Sco X-1 kilohertz QPO frequencies 

Fig. p^ shows the Vj^-p-u-^ relations given by the sequences 
of geodesies around stars constructed using EOS UU and 
EOS C that fit best the v\ and U2 data for Sco X-1. 

The best-fitting sequence for stars constructed using 
the UU EOS is for a nonrotating star with M = 1.89 M©. 
This fit gives x^/dof = 1.4, which is unacceptable at the 
94% confidence level. To achieve fits this good, the sequence 
must include HEGs with ra/rp — 2.1. Such highly eccentric 
geodesies are required in order to generate i^r-i^p^ relations 
that are sufficiently flat. Physically, clumps must be formed 
or injected on a special sequence of highly eccentric geodesies 
with a common periastron radius that is 40% larger than the 
stellar radius and 15% smaller than risco (as shown in Fig. 0, 
if rp > risco is required, as suggested by SV, the fits would 
be much poorer). The best-fitting mass is about 0.3 M© less 
than the maximum stable mass and hence the fit is relatively 
insensitive to the stellar spin rate (x'^/dof would be about 
the same for any spin frequency < 600 Hz). This geodesic 
sequence is the best-fitting sequence for any EOS that can 
support a nonrotating star with M > 1.89 Mq and an equa- 
torial radius Ra < 5.3M. 

The best-fitting geodesic sequence for stars constructed 
using the softer EOS C is for a nonrotating star with 
M = 1.86 Mq, the maximum stable mass for this EOS. The 
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Figure 3. Upper panel: Vr-v^ relations given by two geodesic 
sequences that best fit the v\ and V2 data for Sco X-1. The sohd 
curve shows the v^-u^ relation given by the geodesic sequence that 
fits best overall (M = 1.89 Mq, v^ = 0, rp = 5.3M; x^/dof = 
1.4). This is the best-fitting sequence for any EOS, such as UU, 
that can support a nonrotating star with M > 1.89 Mq and an 
equatorial radius Re < 5.3M. The dotted curve shows the Ur- 
j/p, relation given by the best-fitting sequence for EOS C (M = 
1.86 Mq, V, = 0, rp = 5.1M; x^/dof = 2.2); the fit would be 
worse for a rotating star (see text). To achieve fits this good, both 
sequences must include highly eccentric geodesies. Lower panel: 
Circumferential apastron radii of the highly eccentric geodesies 
that give the frequencies displayed in the upper panel. 



smaller maximum mass of C stars requires geodesies that are 
more eccentric than for UU stars, in order to keep i/r within 
the observed range of Au. This causes the slope of the best- 
fitting t'r-z^K relation to be flatter than for UU stars and flat- 
ter than the observed correlation (x^/dof = 2.2). Requiring 
the star to have a spin rate ^ 300 Hz, hke that expected 
for the neutron stars in the kilohertz QPO sources, would 
increase all orbital frequencies, making the fit significantly 
worse (the maximum mass is not increased significantly for 
the relevant spin rates). We conclude that the geodesic pre- 
cession hypothesis is inconsistent with the kilohertz QPO 
frequencies observed in Sco X-1 if neutron star matter is as 
soft as EOS C. 



I GjY 340 I kilohertz QPO frcqucnoica 



3.2.2 

Fig. y shows the v^-v.^^ relations given by the sequences of 
geodesies for neutron stars constructed using EOS C and 
EOS UU th at best fit the v\ a nd V2 data for the Z source 
GX 340-1-0 djonker et al. 199^ ). The fit is again better for 
higher mass stars, because the observed frequencies of the 
upper kilohertz QPO in GX 340-1-0 are relatively low when 
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Figure 4. Upper panel: Vr-i'^^ relations given by two geodesic 
sequences that best fit the v\ and V2 data for the Z source 
GX 340-1-0. The solid and dotted curves show respectively the u,-- 
v^ relations given by the best-fitting sequences for neutron stars 



p = 7.0M; 
p = 5.8M ; 



constructed using EOS UU (M = 2.12 M©, u^ = 0, 

X^/dof = 0.84) and EOS C (M = 1.86 Mq, 1^^ = 0, 

X^/dof = 2.3). Both sequences require highly eccentric geodesies 

(''a,max/''p = 1.6-2.8). Lower panel: Circumferential apastron 

radii of the highly eccentric geodesies that give the frequencies 

displayed in the upper panel. 



interpreted as azimuthal frequencies near the ISCO. Hence 
stars constructed using EOS C fit significantly worse than 
stars constructed using EOS UU, because the former have 
smaller maximum stable masses. 

The best-fitting geodesic sequence for stars constructed 
using EOS C is for a nonrotating star with a mass equal 
to the maximum stable mass for this EOS, requires a very 
high eccentricity (r-a,max/7'p = 2.8), and does not fit the data 
adequately (x^/dof — 2.3), even though the uncertainties in 
the measured frequency differences /S.v axe rather large. The 
disagreement would be worse for a rotating star. Hence we 
conclude that the geodesic precession hypothesis is inconsis- 
tent with the kilohertz QPO frequency correlation observed 
in GX 340-1-0 if neutron star matter is as soft as EOS C. We 
caution, however, that a better understanding of the uncer- 
tainties in measurements of V2 above 700 Hz caus ed by the 
higher n oise power density at these frequencies (Jonkeret 



al. 1998) could affect this conclusion. 



The best-fitting geodesic sequence for the UU EOS 
again requires highly eccentric (ra,max/''p = 1-6) geodesies 
with a fixed periastron radius more than twice as large as the 
star's equatorial radius. The sequence of geodesies around a 
nonrotating star that fits best has a mass only about 0.1 Mq 
less than the maximum stable mass for this EOS; x^/dof for 
the best-fitting sequence is 0.84. The best-fitting geodesic se- 
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quence and the x^/dof would be approximately the same for 
any spin rate < 200 Hz. 



3.2.3 4U 1 728- 34 kilohertz QPO frequencies 

The observed frequencies of the kilohertz QPOs in the atoll 
source 4U 1728—34 narrow further the range of equations 
of state consistent with the geodesic precession hypothesis. 
Fig. H shows the Vi-v^i relation given by the sequence of 
geodesies for stars constructed using the UU EOS that best 
fit the v\ and V2 data for 4U 1728—34. This fit assumes 
that the stellar spin frequency is the 363 Hz frequency of 
the burst oscillations in this source. The fit would be very 
similar for any spin frequency < 500 Hz. 

Again, a fr-J^K relation qualitatively similar to the 
observed v\-i/2 correlation can be achieved only by us- 
ing HEGs. The I'r-J^K relation given by the best-fitting 
geodesic sequence for stars constructed using the UU EOS 
has ra,max/^p = 2.7 and gives x^/dof — 1.7. The large 
scatter in the observed l\v-U2 relation for upper kilo- 
hertz QPO frequencies in the range 1150 ^ 1^2 ^ 1200 Hz (see 
Mendez & van der Klis 1999) contributes significantly to the 
X^/dof. This scatter appears inconsistent with the simple, 
smooth frequency relations characteristic of geodesic pre- 
cession models. The best-fitting geodesic sequence for stars 
constructed using the modern, realistic EOS A18-|-UIX-|-5ui, 
gives the same model parameters, except that the best- 
fitting value of Tp is 1.05i?c instead of 1.09iZc (see Table H 
below) . 

Stars constructed using harder equations of state, such 
as the early EOS L, which are now thought to be unreal- 
istically hard, have larger maximum masses but also larger 
equatorial radii. This eliminates some of the highly eccen- 
tric geodesies that are available if the EOS is softer, mak- 
ing the geodesic precession hypothesis inconsistent with the 
kilohertz QPO data for such hard equations of state. The 
4U 1728—34 frequency data illustrate this. The flatness of 
the /^v-v.^ correlation for 1^2 < 1100 Hz (see Fig. pi) con- 
strains the mass of the neutron star to be <; 2Mq. The large 
equatorial radii of stars constructed using hard equations 
of state leave room only for less highly eccentric geodesies, 
which produce z^r-i^K relations that fall steeply at high z/j^ 
and are therefore inconsistent with the observed A!/-i/2 cor- 
relation. 



Our findings up to this point show that //^ 



relations 



similar to the i'i-U2 correlations observed in the kilohertz 
QPO sources require highly eccentric geodesies and are not 
possible for geodesies around neutron stars constructed with 
equations of state that are either significantly softer or signif- 
icantly harder than the UU EOS, i.e., the geodesic precession 
hypothesis is inconsistent with softer or harder equations of 
state. Hence, from now on we consider only EOS UU. The 
realistic modern EOS A18-|-UIX-|-^i;i, gives results that differ 
only very slightly from those obtained with the UU EOS. 

3.2.4 GX 5-1 and GX 17+2 kilohertz QPO frequencies 

Fig. o shows the ten pairs of kilohertz QPO frequencies mea - 
sured so far in the Z source GX 5—1 ( Wijnands et al. 1998 ). 
Most have relatively large uncertainties and one (the cir- 
cled point in the figure) appears discordant. If all ten fre- 
quency pairs are included, the v^-v^ relation given by the 
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Figure 5. Upper panel: Vi-^y. relations given by two geodesic 
sequences that best fit the V1-V2 data for the atoll source 
4U 1728—34. The solid and dotted curves show respectively the 
fc'r-t'j^ relations given by the best-fitting sequences for neutron 
stars constructed using EOS UU (M = 1.76 M©, r^ = 4.7M = 
1.09Ra = 12.4km; x^/dof = 1-7) and EOS L (M = 1.97 Mq, 
rp = Re = 5.3M = 15.4km; x^/dof = 8.6), assuming u^ = 
363 Hz, the burst oscillation frequency. Both sequences require 
highly eccentric geodesies (ra,max/fp = 2.7 for the best-fitting 
UU sequence). Lower panel: Circumferential apastron radii of the 
highly eccentric geodesies that give the frequencies displayed in 
the upper panel. 



best-fitting sequence of geodesies for stars constructed using 
EOS UU has a rather large x'^/dof; discarding the discor- 
dant frequency pair gives a smaller x'^/dof of 1.2. However, 
values of x^/dof this small are possible only if the star has a 
mass M ~ 2.1-2.2 Mq near the maximum stable mass and 
highly eccentric geodesies are included (ra,max/''p ~ 2.1). 

Only eight pairs of kilohertz QPO frequencies have been 
measured so far in the Z source GX 17-1-2 (Wijnands et al. 
1997; see Fig. hll below). The fr-fj^ relation given by the 
best-fitting sequence of geodesies for stars constructed using 
EOS UU is acceptable (x^/dof ^ 0.4), but only if all the 
geodesies are highly eccentric (ra.max/rp ~ 3-4). However, 
the uncertainties in these frequency measurements are large. 



3.3 Other geodesic frequencies 

We have also explored whether the frequencies of the two 
kilohertz QPOs are consistent with geodesic frequencies 
other than i/^p and v^^. As discussed in Section m the 
most prominent frequencies generated by clumps moving 
on geodesies around neutron stars are j/r, 21^1, v^ + Vr, and 
i^j^ -|-2!/r. We find that for realistic neutron star models, these 
frequencies are not in the same range as the frequencies of 
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Figure 6. Upper panel: Vi-v^ relations given by two geodesic 
sequences that best fit the v\-V2 data for the Z source GX 5—1. 
The solid and dotted curves show respectively the v^-v^ relations 
given by the best-fitting sequences for neutron stars constructed 
using the UU EOS if all the frequency measurements are included 
(M = 2.17 M0, rp = 6.5M; x^/dof = 3.6) or, alternatively, if 
all but the circled measurement are included (M = 2.09 Mq, 
Tp = 6.5M; x'^/dof = 1.2). Fits this good require highly eccentric 
geodesies with ra/rp fs 2. Lower panel: Circumferential apastron 
radii of the highly eccentric geodesies that give the frequencies 
displayed in the upper panel. 



the kilohertz QPOs. The two frequencies v^ and i/-^ + 1/^ are 
in the same range as the frequencies of the kilohertz QPOs 
only for unrealistically hard neutron-star equations of state, 
such as EOS L; however, even for such hard equations of 
state, the best-fitting geodesic sequences give frequency re- 
lations that have qualitatively different shapes from those 
observed and give xVdof > 10 for GX 340+0 and GX 5-1. 
We therefore conclude that geodesic frequencies other than 
Vp^-p and v-^ are not viable explanations of the frequencies of 
the kilohertz QPOs. 



3.4 LFQPO and HBO frequencies 

In the previous two sections we investigated whether se- 
quences of geodesies can be found that give J'ap'^k relations 
consistent with the kilohertz QPO frequencies observed in 
the Z and atoll sources. There we showed that the best- 
fitting sequences of infinitesimally or moderately eccentric 
geodesies give J^^p'^k relations qualitatively different from 
the kilohertz QPO frequency correlations in all the sources 



to the kilohertz QPO frequencies observed in the sources 
considered are possible only if highly eccentric geodesies 
{fa/rp ~ 2-4) are included in all the sequences. HEG se- 



consistent with the i'i-U2 correlations observed in Sco X-1 
and 4U 1728—34, although the values of x'^/dof are not for- 
mally acceptable. HEG sequences can also be found that 
formally fit the fi-u^ correlations observed in GX 17-1-2 and 
GX 340-1-0, but there are only a few frequency measure- 
ments for each these sources and the uncertainties are large. 
However, the geodesic sequence that best fits the GX 5—1 
frequency data appears to be inconsistent with these data. 

SV and SVM proposed that gas orbits the neutron stars 
in the Z and atoll sources in the form of clumps moving on 
slightly inclined geodesies and that the low-frequency peaks 
observed in power spectra of the Z and atoll sources are 
caused by nodal precession of these geodesies. As explained 
in Section hi, i/j^p is too low to explain the frequencies of 
these low-frequency features. SV and SVM therefore sug- 
gested that clump nodal precession somehow produces an X- 
ray oscillation with frequency 2!/j^p . (However, as explained 
in Section H, we were unable to find a clump mechanism that 
generates a significant peak at 2i/pjp.) 

Here we follow up the proposal that the low-frequency 
peaks are generated by nodal precession by investigating 
whether sequences of geodesies can be found that give J^^p" 
j/^ and 2i/^p-iy^ relations that are both simultaneously con- 
sistent with the i^L'^^z and i^i-i/2 correlations observed in 
the Z and atoll sources. In contrast to the (^^p'^'k relations 
studied in the previous section, the 2i'^p-i/^ relations stud- 
ied here are sensitive to the stellar spin rate, because of the 
gravitomagnetic effect. 

Before proceeding, we caution that measurements of the 
centroid frequencies of the bumps observed in atoll-source 
power spectra can be affected significantly by systematic er- 
rors that arise from uncertainties in modeling the noise con- 
tirmum in the vicinity of these features (M. Mendez and M. 
van der Klis 1999, private communication). The systematic 
errors are often difficult to estimate and typically are not 
reported in the literature. Instead, only statistical errors are 
reported. We use the reported errors to make a preliminary 
assessment of whether the fits below are acceptable, but one 
should bear in mind that the total errors may be somewhat 
larger. 



3.4.1 4U 1728-34 LFQPO and kilohertz QPO frequencies 

The 2i/jjp-i^j^ relation given by the geodesic sequence which 
gives the f^p'^K relation that best fits the v\-i'2 data in 
4U 1728—34 (see Section |!^ and Fig. H) has a shape that is 
qualitatively different from the observed Vi^-v^ correlation 
in this source (see Fig. m. The predicted relation shown 
is for a stellar model constructed using the UU EOS and a 
spin rate equal the 363 Hz burst oscillation frequency. The fit 
cannot be improved significantly by increasing or decreasing 
the spin rate. 

We also considered stellar models constructed using the 
modern realistic EOS A18-|-UIX-|-5?;6 (see Section ^^ . The 
best-fitting stellar mass is 0.1% larger and the dimensionless 
angular momentum is 3% larger {j = 0.156 vs. j — 0.149) 
for this EOS, because it is slightly harder. This increases 
the nodal precession frequency for geodesies with a given 



i^^, but only by ^4%. 

The predicted 2i/j^p-i/j^ relation shown in Fig. tn is clearly 
inconsistent with the observed J^l-^'2 correlation: the pre- 
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Figure 7. Upper panel: Comparison of tlie ui^-uk behaviour 
seen in 4U 1728—34 with the 2ut^p-ui,^ relation predicted by the 
geodesic sequence that gives the best fit to the i'i-i'2 data for this 
source (see Fig. 5). Lower panel: Circumferential apastron radii of 
the highly eccentric geodesies that give the frequencies displayed 
in the upper panel. 



dieted relation does not reproduce the abrupt drop in v^ 
by a factor > 2 that occurs at ^2 ~ 900 Hz in the observed 
correlation, nor is it consistent with the observed z^^'^^s '-or- 
relation either before or after the abrupt drop. The predicted 
^NP"^K relation misses many of the observed values of u^ by 
a factor ~2. This qualitative disagreement between the fre- 
quency behavior predicted by the nodal precession hypoth- 
esis and the frequency behaviour observed in 4U 1728—34 
occurs for any realistic neutron star model and spin rate. 

Very recent analyses of 4U 1728—34 power spectra (Di 
Salvo et al. 2000, private communication) indicate that the 
break frequency of the low-frequency broad-band noise com- 
ponent and the frequency of the LFQPO are related: the 
frequencies of the break and the LFQPO both increase as 
the frequency of the upper kilohertz QPO increases from 
400 to 900 Hz, but at this point the existing LFQPO dis- 
appears and the noise feature peaks up into a new LFQPO. 
The frequency of the old and new LFQPOs are not harmon- 
ically related. This behaviour appears inconsistent with any 
simple interpretation in terms of nodal precession of HEGs. 



3.4.2 GX 5—1 HBO and kilohertz QPO frequencies 

The behaviour of 2i/^-p is also qualitatively inconsistent with 
the behaviour of the HBO frequency in GX 5—1. This is 
illustrated in Fig. fel which compares the !/hbo-'^2 correlation 
observed in GX 5—1 with the 2v„-, 



400 



N 



^ 300 



200 
60 

^ 50 

^ 40 

^ 30 






20 

10 



20 



^ 10 



- 


1 — 1 


^ " 




- 


^ y^ ■ 




H — 


^<N ' 




V \\: 


- , , , \1 


_ 1 1 1 i 


y 


y 


— y 


/ 


y 


, ^ 


^ 


y 


~ ^ ^^ 




^ . ' ^ — 








^ ^ ^ ---' — ^ 


^-^'' 


^_ — - 


1 1 1 1 


1 1 1 T 


$^ 


N-^- 


X^"- 


^"^ 


^Sy^^"*!^^^^ 






^^^^^^- 


^ 1 1 1 ^ 



400 



600 

V 



K' 



v.. 



800 

[Hz] 



1000 



relations predicted 



Figure 8. Geodesic sequences around UU stars that best fit the 
1^1-1^2 data for GX 5—1, for different stellar spin rates. The solid, 
long-dashed, and short-dashed curves show the relations given 
by the geodesic sequences that fit best for v^ = 100 Hz (M = 
2.20 Mq, rp = 6.3M; x^/dof = 3.6), u^ = 363 Hz (M = 2.21 Mq, 
Tp = 5.8M; x^/dof = 4.5) and u^ = 600 Hz (M = 2.22 Mq, 
rp = 5.3M; x^/dof = 5.6). Re = 3.02M for all three stellar mod- 
els. Top panel: Comparison of the best-fitting ^^ap'^^k relations 
with the observed t'i-!^2 correlation. Middle panel: Comparison of 
the predicted i^^p'^^k relations with the observed t'HBO-i'2 corre- 
lation. Bottom panel: Circumferential apastron radii of the highly 
eccentric geodesies that give the frequencies plotted in the upper 
and middle panels. 



by several geodesic sequences constructed to give J^ap'^k 
relations that fit as well as possible the 1^1-1^2 data. 

The basic difficulty is that in order to construct a se- 
quence of HEGs that gives a J^^p'^^k relation similar to the 
observed !/i-i^2 correlation, the stellar spin rate must be 
^400 Hz, but for spin rates this low, 2u^p is much smaller 
than the observed HBO frequencies. If the spin rate is in- 
creased to the values > 600 Hz required to lift 2i/^p into the 
observed frequency range of the HBO, i^j^ is shifted upward 
so much that the J'ap'^k relation given by the best-fitting 
geodesic sequence is much flatter than the observed vuBO-t^K 
correlation. Furthermore, attempting to fit this correlation 
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Figure 9. Geodesic sequences around UU stars that best fit the 
v\-V2 data for Sco X-1, for different stellar spin rates. The solid, 
long-dashed, and short-dashed curves show the relations given 
by the geodesic sequences that fit best for v^ = 450 Hz (M = 
2.13Mq, Re = 3.32M, rp = 4.9M), u^ = 500 Hz (M = 2.16 Mq, 
iJo = 3.27M, Tp = 4.8M) and u^ = 550 Hz (M = 2.18 Mq, 
Re = 3.19Af , rp = 4.8M). x^/dof = 1.4 for all three sequences. 
Top panel: Comparison of the best- fitting J^ap'^^k relations with 
the observed i'i-U2 correlation. Middle panel: Comparison of the 
predicted 2f„p-Vj^ relations with the observed i^bo-i^2 correla- 
tion. Bottom panel: Circumferential apastron radii of the highly 
eccentric geodesies that give the frequencies plotted in the upper 
and middle panels. 



drives the stellar mass to the maximum stable mass for the 
spin rate considered, which is implausible. 



3.4-3 Sco X-1 HBO and kilohertz QPO frequencies 

The behaviour of 2^^^ is qualitatively inconsistent with the 
observed behaviour of the HBO frequency in Sco X-1. This 
is demonstrated in Fig. ^ which compares the iyiiBO-1^2 cor- 
relation observed in Sco X-1 with the 2!/„p-i^j^ relations pre- 
dicted by three sequences of HEGs constructed to give the 
!/.„-;/„ relation that best fits the 1/1-1/2 data for Sco X-1, for 
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Figure 10. The geodesic sequence around UU stars that best 
fits the 1/1-1/2 data for GX 17-1-2. The best-fit parameters are 
M = 1.82 Mq and rp = 4.2M fa 1.01 /?„; X^/dof = 0.39. Top 
panel: Comparison of the best-fitting i/r-t/j^ relation with the ob- 
served A1/-U2 correlation. Middle panel: Comparison of the pre- 
dicted 2Uj^.p-u^ relation with the observed vb.bo-i/2 correlation, 
for v^ = 550 Hz. Bottom panel: Circumferential apastron radii of 
the highly eccentric geodesies that give the frequencies plotted in 
the upper and middle panels. 



various assumed spin rates. For the reason discussed previ- 



ously, the predicted i/j^ 



relations are insensitive to the 



star's spin frequency for the frequencies considered and they 
therefore lie almost on top of one another. 

Here the basic difficulty is that the 2i/j^p-i/j^ relations 
predicted by the geodesic sequences that give acceptable fits 
to the i/\-i/2 data rise much more steeply than the observed 
i/HBO-1^2 correlation, regardless of the spin rate of the star. 



3.4.4 GX 17+2 HBO and kilohertz QPO frequencies 

The behaviour of 2i/^p predicted by the sequence of HEGs 
that best fits the 1/1-1/2 data for GX 17-f2 is qualitatively 
similar to the i/ubo-'-'2 correlation observed in this source, 
but is inconsistent in detail. This is shown in Fig. hd. The top 
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Figure 11. The geodesic sequence around UU stars that pro- 
vides the best simultaneous fit to the wi-U2 and j^hbO'^^k data 
for GX 17+2. The best-fit parameters are M = 2.21 M0, Vp = 
5.0 A/ « 1.7iic, and u^ = 660 Hz; x^/dof = 4.3. The open cir- 
cles show the frequency triplets {f^p ; 1 !^k i ' ^'^np i } that fit best. 
Top panel: Comparison of the resulting Vr-v^^ relation with the 
observed Av-V2 correlation. Middle panel: Comparison of the re- 
sulting 2fjjp-j/j, relation with the observed vhbo-^2 correlation. 
Bottom panel: Circumferential apastron radii of the highly eccen- 
tric geodesies that give the frequencies plotted in the upper and 
middle panels. 



panel compares the Vr-v^. relation given by the best-fitting 
sequence of HEGs with the observed A:/-!^2 correlation. Very 
highly eccentric geodesies (^a/rp « 4) are required in order 
to obtain an acceptable fit. The x^/dof is substantially less 
than 1 because of the large uncertainties reported for the u\ 
and V2 measurements. The fit is relatively insensitive to the 
spin frequency, for low to moderate spin frequencies. 

The middle panel of Fig. hd compares the 2v-^-p-v-^ re- 
lation given by the sequence of HEGs used in the top panel 
that is most similar to the observed !^bo-'^2 correlation. The 
relation shown is not a joint fit to the HBO and kilohertz 
QPO frequencies, but was instead obtained by adjusting the 
spin frequency to give the 2Vj^p-v^ relation that is closest to 
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Figure 12. Upper panel: Comparison of the 2v-^-p-Vj^-p rela- 
tions given by sequences of infinitesimally eccentric geodesies with 
the correlation between the frequencies v^ and v^ observed in 
Cir X-1. The solid and dotted curves show respectively the fre- 
quency relations given by sequences of MEGs around UU stars 
that give relations similar to the v^^-v^ correlation observed at 
u^ < 100 Hz, for Us = 700 Hz (M = 2.1 Mq, j = 0.26) and 
Us = 900 Hz (M = 1.3 Mq, j = 0.50). The dashed and dash- 



dotted curves show respectively the 2u^ 



relations given by 



moderately eccentric geodesies around black holes with the same 
values of M and j as for the high- and low-mass neutron stars. 
Lower panel: Circumferential radii of the neutron star geodesies 
that give the frequencies plotted in the upper panel. 



the observed V11BO-V2 correlation. The spin frequency that 



gives the best fit is 550 Hz. Although the 2v^ 



relation 



passes near the j^hbo-J^2 data, it misses many of the data 
points by much more than their uncertainties. The best- 
fitting Vj^p-V2 relation clearly does not correctly predict the 
t^HBO-'^2 correlation observed. 

We then explored whether an acceptable joint fit to the 
1^1-1^2 and VUS0-V2 correlations is possible. Fig. hll shows 
the 2u^-p-v-^ and i'ap'^k relations given by the geodesic se- 
quence that best fits both the t'HBO-i'2 and the U1-U2 data. 
The high precision of the HBO frequency measurements 
forces the fit to make the 2i/jjp values agree closely with 
them. However, this is possible only if the eccentricity of 
the geodesies is reduced, but this makes the J^ap'^^k rela- 
tion steeper than the observed v\-V2 correlation. The best-fit 
mass is driven to the maximum stable mass in an effort to 
compensate, but the compensation is insufficient to give an 
acceptable fit. Also, the periastron radius is driven far from 
the stellar surface (rp = 1.7 -Re). Although the geodesies 
favoured by the attempt to fit jointly the !^bo-J^2 and Vi_-U2 
data are less eccentric than those required by the V\-V2 data, 
they are still highly eccentric [r^/vp ~ 2). 
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Figure 13. EOS UU-based neutron star models (open circles) 
tiiat allow 2fjjp-i/^p relations similar to the v^^-v^ correlation 
observed in Cir X-1 for v^ < 100 Hz. The solid line in the upper 
panel shows for comparison the j-M track along which sequences 
of geodesies around a spinning black hole can be constructed that 
give 2Uj^-p-Vp^-p relations similar to the observed v^-v^ correlation. 



3.5 Cir X-1 QPO frequencies 

Although the peculiar accreting X-ray star Cir X-1 is gen- 
erally thought to be a neutron star (see Tennant, Fabian 
& Shafer 1986), Psaltis et al. (1999b) have suggested that 
its QPOs may be related to some of the QPOs observed in 
black-hole candidates, because power spectra of its X-ray 
emission are similar in some ways to the power spectra of 
black holes. SVM have proposed that the frequencies u-^ and 
y^ of the higher- and lower-frequency features observed in 
Cir X-1 are the apsidal precession frequency j/^p and twice 
the nodal precession frequency v^^-p of a sequence of slightly 
inclined lEGs around a neutron star, although they did not 
attempt a fit of the resulting frequency relations with the 
frequencies of the QPOs observed in Cir X-1. 

Here we first investigate whether a sequence of lEGs 
can be found that gives a 2v-^-p-Vj^^p relation that is consis- 
tent both with the QPO frequencies reported by Shirey et 
al. (1996, 1998) and Shirey (1998), and with the centroid 
frequencies of the peaked noise and QPO reported by Ten- 
nant (1987). No error bars were reported by Shirey et al., 
so our conclusions are necessarily tentative. The solid and 



dashed lines in Fig. 12 show the 2v-^ 



relations given by 



the two sequences of lEGs that roughly bracket the range 
of masses and spin frequencies that agree qualitatively with 
the observed v-^-v-^ correlation for v-^^ < 100 Hz (see Fig. FL3). 
In order to achieve qualitative agreement with the observed 
frequency correlation, the neutron star must be rapidly spin- 
ning (u^ > 650 Hz) and the radii of the geodesies in the se- 
quence must vary smoothly by a factor of 2.5, from ~ 25 M 
to ~ 10 M, as the frequencies increase. 

Inclusion of the higher-frequency {v-^ ~ 12 Hz, v-^ ~ 
200 Hz) measurement reported by Tennant (1987) makes it 
much more difficult to find a sequence of geodesies that 




AP 



[Hz] 



relations 



Figure 14. Upper panel: Comparison of the "iv^^ 
given by sequences of highly eccentric geodesies with the correla- 
tion between the frequencies v^^ and v^ observed in Cir X-1. The 
solid and dotted curves show respectively the frequency relations 
given by sequences of HEGs around UU stars that give relations 
similar to the ^i^-v^ correlation observed at v^ < 100 Hz, for 
Us = 550 Hz (M = 2.2 Mq, j = 0.18, rp = 4.1Af) and Us = 900 Hz 
(A/ = 1.3 Mq, j = 0.50, rp = 11. OM). The dashed and dash- 



dotted curves show, for comparison, the 2u^ 



relations given 



by geodesies around black holes with the same values of M, j, 
and Tp as for the high- and low-mass neutron stars. Lower panel: 
Circumferential radii of the neutron star geodesies that give the 
frequencies plotted in the upper panel. 



gives a frequency relation similar to the observed correla- 
tion. Only lECs around stars with very low masses and very 
high spin rates give frequency relations that pass close to 
the lower-frequency data of Shirey et al. (1996, 1998) and 
then bend over sufficiently to pass near the higher-frequency 
measurement. Sufficient flattening of the 2i/j^p-i^^p relation 
is achieved only for stars that are highly deformed by their 
rapid spins and therefore have large mass quadrupole mo- 
ments whose effect largely cancels (at small radii) the pro- 
grade precession produced by frame-dragging. The effect of 
a large mass quadrupole may be seen clearly in Fig. [12 by 
comparing the 2i^j^p-i^^p relations for neutron stars with 
the relations for black holes with the same M and j; no- 
tice the particularly large difference for M = 1.3 Mq and 
j — 0.50; the neutron star model is still fairly far from the 
mass-shedding limit, which is 0.93 Mq. 

Given the difficulty of approximating the J^l'^^h correla- 
tion observed in Cir X-1 using lEG sequences, we explored 
whether it would be easier to reproduce the correlation us- 
ing a sequence of HEGs. Our results are shown in Fig. |l4[ 
The lower- frequency portion of the J^l"^h correlation drives 
inclusion of very highly eccentric [r^jr^ > 6) geodesies that 
extend far from the neutron star (ra/-Ro > 10) in the se- 
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Table 3. Fits of geodesic sequences to the measured frequencies of the kilohertz and lower-frequency QPOs. 
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UU 
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3.02 


5.62 
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HEG 


UU 


600 


2.22 
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9.9 


3.02 


5.37 


5.3 


7.6 
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5.6 
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UU 
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12.4 
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"Number of frequency measurements (n) minus the number of parameters in the model (p). ''Fits for nonrotating stars do not depend on the EOS, 
provided that it is hard enought to support the best-fit stellar mass and that R^ < fp,min- The values of R^ listed for all nonrotating stars are for 
models constructed using the UU EOS, except for the 1.86 M© stars, which were constructed using EOS C. '^These models are not fits and would 
give very large values of x^/dof. They are the geodesic sequences reported by SV as adequately fitting the Sco X-1 kilohertz QPO frequency data 
(see text). ''No innermost stable circular orbit exists, because the radius of the star is too large. "^No formal fits were possible, because no errors 
have been reported for most of the frequency measurements used; the fit was therefore done by eye, giving each point the same weight. 



quence of HEGs. Inclusion of the higher- frequency measure- 
ment of Tenant (1987) again makes it much more difficult 
to find a sequence of geodesies that gives a frequency rela- 
tion similar to the observed correlation. Only sequences of 
HEGs around stars with very low masses and very high spin 
rates give frequency relations that pass near all the mea- 
surements. The geodesic sequences that work best have in- 
finitesimally or moderately eccentric geodesies as their most 
compact members, but even these geodesies have periastron 
radii much larger than the large equatorial radii of these 
stellar models {vp/Rc = 1.7). 



4 POWER DENSITY SPECTRA 

SV, SVM, and Karas (1999) simply assumed that the most 
prominent frequencies in the X-ray waveforms produced by 
clumps moving on geodesies around a neutron star are 2i'j^p , 
r/^p , and f^ . No attempt was made to compute the frequen- 
cies that would be generated or the relative amplitudes of the 



oscillations at these frequencies. Here we report calculations 
of the frequencies and amplitudes that would be produced. 
In order to determine which are the dominant frequen- 
cies that would be generated and the relative amplitudes of 
the oscillations at these frequencies, we carried out an exten- 
sive suite of time-dependent, numerical simulations of the 
geodesic motion of clumps around spinning neutron stars. 
For the purposes of this investigation we simply assumed 
that all clumps are formed and move on the same geodesic, 
that they are held together by some force, and that they 
move forever on this same geodesic. We modeled the clumps 
as spheres of constant radius d. In order to facilitate the com- 
putations, we assumed d = 0.17?e, which is about 100 times 
larger than allowed by the dynamical constraints (see Sec- 
tion n). We treated the motion of the clumps in full general 
relativity and took into account special relativistic Doppler 
and other effects on the radiation absorbed and emitted by 
the clumps, but we did not perform ray tracing. Results 
using full ray tracing would be qualitatively similar to the 
results reported here. We surveyed a broad range of geodesic 



© 2000 RAS, MNRAS 000, 000-000 



Eccentric orbits and QPOs in neutron star X-ray binaries 17 



0.01 



0.001 



0.01 



^ 0.001 



0.0001 
0.1 



K 001 



0.001 



0.1 



0.01 



0.1 - 



0.01 



0.01 



0.001 



a 



H \ H 



d 



l/„+2!/. 



2u,. 



H \ k H h 



1 1 ''" I 1 h 



H ^ 



H 1 ' " l|l 



H ^ 



' '' I' I ^ ^ 



-H h 



200 



400 



600 



800 



1000 
1/ 



1200 



1400 



1600 



1800 



Figure 15. Typical power spectra produced by six different effects, for rouglily spherical clumps on highly eccentric, tilted geodesies 
around a spinning neutron star. The heights of the peaks are not proportional to their total power, because their widths are different 
(total powers are listed in Table H. From top to bottom, the effects considered are (see text for details): (a) occultation of a uniformly 
bright neutron star surface; (b) occultation of a belt with a Gaussian brightness distribution of semi-width O.lRo around the stellar 
equator; (c) reflection of radiation from a uniformly bright neutron star surface; (d) radiation by a clump with a radius-independent 
proper luminosity; (e) radiation by a clump with a radius-dependent proper luminosity; (f) radiation by a clump in a tilted orbit that 
emits only while it is passing through a rarified disc of height h{r) = 0.05 r in the equatorial plane of the system. In this example, 



M = 1.76 Mq, u^ = 363 Hz rp = 4.7Af, and r-a = 10. OM, yielding u^ = 1107Hz, Vr = 343 Hz, and u^ 



19 Hz. The geodesic was 



assumed to be inclined 10° relative to the equatorial plane. The line of sight was assumed to be inclined 45° . As explainlained in the 
text, the actual power spectrum produced by a particular clump model is a linear combination of two or more of the spectra shown. 



parameters, including rp, ra, and the tilt of the geodesic, and 
a wide range of viewing angles. 



4.1 Highly eccentric geodesies 

We simulated six different effects that could modulate the X- 
ray flux seen by a distant observer. The six panels of Fig. M 
show power spectra that are typical of the spectra produced 
by each of these effects, for clumps moving on geodesies 
with the high eccentricities required to produce kilohertz 
QPO frequencies qualitatively similar to those observed in 



the atoll and Z sources. The heights of the peaks are not pro- 
portional to their total power, because the motions involved 
are not periodic, which causes the widths of the peaks to 
differ. The frequencies and relative total powers of the most 
prominent peaks are listed below in Table kl As discussed 
further below, the actual power spectrum produced by a par- 
ticular clump model is a linear combination of two or more 
of the spectra shown in Fig. hs. We now discuss briefly each 
of six effects (the labels a-f used here match the labels of 
the panels in Fig. hsl). 

(a) Occultation of the stellar surface. An optically thick 
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Figure 16. Typical power spectra produced by the same six effects considered in Fig. 15, but for moderately eccentric geodesies. As 
in Fig. 15, the heights of the peaks are not proportional to their total power, because their widths are different (total powers are listed 
below, in Table M). The stellar mass and spin rate and the inclinations of the geodesic and line of sight are the same as in Fig. 15, but 



5.6M and ra = 6.7M, corresponding to an eccentricity e = 0.09 and yielding u^^ = 1211Hz, Ur = 358 Hz, and h'j^ 



19 Hz. As 



explainlained in the text, the actual power spectrum produced by a particular clump model is a linear combination of two or more of 
the spectra shown. 



clump moving around a neutron star on a geodesic may 
occult a portion of the stellar surface, thereby modulat- 
ing the X-ray flux from the surface seen by a distant ob- 
server. Occultation occurs only for geodesies that satisfy 
Tp < {Ra + d)/ cosi, i.e., that closely approach the stellar 
surface. Some of the best-fit geodesies listed in Table do 
not satisfy this constraint unless the observer's sightline is 
close to the plane of the disc. 

In the absence of apsidal and nodal precession, the star 
would be occulted periodically, with frequency u^^ — fr. 
However, general relativistic effects near the neutron star 
make v^ > v,.- For the geodesies considered here, a clump 
typically orbits the star three or four times each time it 
moves radially inward and outward at the radial epicyclic 



frequency v^ . The rate of azimuthal phase advance varies sig- 
nificantly with the radial position of a clump and hence the 
occultation caused by the azimuthal motion is not periodic. 
The only strictly periodic motion is the radial epicyclic mo- 
tion, which causes the observed fiux to vary with frequency 
Ur as the extent of the occultation increases and decreases 
with the motion of the clump toward and away from the 
star (depending on the value of r^ and the inclination of the 
geodesic and the line of sight, a clump that is near ra may 
not occult the star at all). 

The strongest peaks in the power density spectrum of 
the X-ray waveform are at the radial epicyclic frequency Ui 
and the combination frequencies i/j^ -I- v^ and f^ + 2ur (see 
Fig. b_aa). The next strongest peak is at 2v-r. Weaker nodal 
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precession sidebands appear at v^+Vr- 
j/^p in this power spectrum and the others discussed below 
only if the geodesic is significantly tilted. There is much less 
power at the azimuthal frequency v-^ and almost none at the 
apsidal precession frequency //^p = Uj^ — Ur- These are the 
frequencies that SV, SVM, and Karas (1999) assumed would 
be the dominant frequencies in the power spectrum. The 
power in the peak at i/^ is typically ~4-5 times less than the 
power in the peak at Vr; t he p ower at Uj^p is too low for the 
peak to be visible in Fig. |l5p . Nodal precession modulates 
the angle between the orbital plane and the observer's line 
of sight, producing a very weak peak at v^^p', there is no 
significant power at 21/^^ . 

(b) Occultation of a bright equatorial band. If radiation 
comes predominantly from a narrow belt around the star's 
rotation equator, nodal precession has a stronger effect on 
the waveform. Fig. [13b shows a typical power spectrum of 
the waveform for this situation. In this example, the surface 
brightness of the star was assumed to vary as exp(— z^/cr^), 
where z is the distance of the point on the stellar surface 
from the equatorial plane and a — O.lRc- 

The strongest peaks are again at Vr, u^ + Ur, and 
Uj^ + 2vi. The next strongest peaks, which are about half as 
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p and v-^ + 2vi — 



strong, are at v-^ - 



and i/„ +2i/r 



. There is a peak 



at Vj^-p that is slightly weaker than these nodal-precession 
sidebands. Again, there is much less power at v^ and al- 
most none at Vj^p. The peak at u-^ is ~3-4 times weaker 
than the three strongest peaks in the power spectrum and 
is significantly weaker even than the peak at v^^p . The next 
strongest peaks, which have almost as much power as the 
peak at v-^^ , are at the nodal sideband frequencies u^ — v^p , 
Vy + i^Np , and v^ + 2fr + v^j. ■ The peak at i^^p is too weak 
to be visible and there is no significant power at 2i/jjp . 

(c) Reflection of stellar radiation by clumps. Fig. FTSc 
shows the power spectrum for a typical waveform generated 
by clumps that reflect radiation coming from the stellar sur- 
face. In this example, we assumed that the surface brightness 
of a clump varies with its distance r from the center of the 
star as 1/r^ and that the flux is proportional to its surface 
brightness and the projected area of the illuminated part of 
the surface seen by the distant observer. The Doppler shifts 
caused by the velocity of the clump relative to the stellar 
surface and to the observer were included. 

The strongest peak produced by refiection is at u-^ +1^1. 
There are also strong peaks at fr and i^^ + 2t'r. There are 
peaks at v^ and Uj^p , but the power in these peaks is ~4-5 
times less than in the peak at u^ + Ur and ~3-4 times less 
than in the peaks at Vy: and v-^^ + 2^^. There is no significant 
power at either z/„p or 2u-^p. 

(d) Radiation by constant temperature clumps. Signifi- 
cant radiation by clumps with the small size and high den- 
sity required by gas dynamical constraints (see Section 0) 
is ruled out by their very small areas compared to that of 
the accreting star and by the upper bound on their tem- 
perature imposed by their origin in a cool disc fiow and the 
requirement that they be in pressure equilibrium with their 
environment (otherwise they will expand and dissipate on a 
timescale much shorter than the lower bound on the dura- 
tions of their wavetrains imposed by the observed coherence 
of the kilohertz QPOs). Nevertheless, for completeness we 
explored the waveform that would be produced by such ra- 
diation. 



Even if the luminosity of a clump is constant in time 
in its rest frame, the general relativistic redshift and other 
effects will cause the luminosity seen by a distant observer 
to increase and decrease with the clump's distance r from 
the center of the star, which varies at the radial epicyclic 
frequency. For the illustrative calculations reported here, 
we took into account only the Doppler shift caused by the 
clump's motion relative to the distant observer. 

A typical power spectrum is shown in Fig. |l5|d. The 
strongest peak is at u-^+u^. There is also a peak at u-^ + 2ut: 
which is about half as strong. There is a peak at j/^p , which 
is also about half as strong as the peak at v-^ -\- Vr, as well 
as a peak at v-^, which is ~4 times weaker than the peak 
at Vy^ + Ur and only about half as strong as the peak at 
u-^ -K 2!^r. There are no other significant peaks. Hence, even 
in this case the peaks at v^ + v^ and v^ + 2vr are as strong 
or stronger than the peaks at v^ and Vj^p . 

(e) Radiation by clumps heated by radiation from or 
interaction with the star. Again, significant radiation by 
clumps is ruled out by their very small areas and by the 
upper bound on their temperature imposed by the require- 
ment that they be in pressure equilibrium. Furthermore, any 
interaction of the clumps with the stellar surface or with a 
stellar magnetic field that is strong enough to release an ap- 
preciable fraction of their kinetic energy would disrupt them 
and cause their orbits to decay. Nevertheless, for complete- 
ness we also explored the waveform that would be produced 
by such radiation. 

The rest-frame luminosity of a clump heated by radi- 
ation from the star or by interaction with a stellar atmo- 
sphere or magnetic field necessarily varies strongly with r. 
The luminosity seen by a distant observer is also affected 
by the variation of the Doppler shift with r. For the illus- 
trative calculations reported here, we assumed that heat- 
ing causes the brightness of a clump to vary with radius as 
l/[{r-R,f + O.OlR',]. 

A typical power spectrum is shown in Fig. fLSe. The 
strongest peak is at v^, with peaks at v^^ + u^ and v-^.^ + 2ur 
that are about half as strong and a peak at 2vy: that is about 
a third as strong. There is a peak at f^, but it is ~ 5 times 
weaker than the main peak and less than half as strong as 
the peaks at v^^ + v^ and v^ + 2ui; there is no significant 
power at u^p , j/^p , or 2u^p . 

(f) Interaction of clumps with a rarefied gaseous disc. 
Finally, we consider the possibility that the luminosity in- 
creases when clumps on tilted geodesies pass through the 
equatorial plane. Such flaring could be caused, for exam- 
ple, by interaction of the clumps with a rarified, remnant 
gaseous disc in the plane near the star. The additional lu- 
minosity produced by enhanced dissipation would have to 
be very small in order to avoid causing the clumps' orbital 
motion to decay more rapidly than is allowed by the ob- 
served coherence of the kilohertz QPOs. This scenario was 
constructed specifically to produce a peak as strong as pos- 
sible at 2//jjp . The results presented here assume the clump 
brightness varies as exp(— z^//i^), where h{r) = 0.05 r is the 
height of the disc, but they are insensitive to the precise 
value of h. 

Even in this scenario, the two strongest peaks are at 
^K + ^r + 2!^„p and v-^ + '^'^n-p ■ The next strongest peaks 
are at v-^ -\- v-^ and v.^+Vi — 2i^pjp and have ~2.5 times less 
power than the peak at i^^ + ^'r + 2i^NP • There are two peaks 
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Table 4. Total power in the prominent peaks in the twelve typical power spectra shown in Figs. 15 and 16.^ 





HEG (rp = 4.7, r^ = 10.0)3 




MEG (rp = 5.6, r^ = 6.7)3 






Frequency ^ 


a b c d e 


f 


abed 


e 


f 


'^NP 


0.10 0.40 — — — 


— 


0.54 0.82 — — 


— 


— 


Si'np 







— 0.24 — — 









Vr + I'nf 
2Vr 



0.90 0.90 0.79 — 1 — 

— 0.21 — — — — 

— 0.23 — — — — 
0.16 0.22 — — 0.30 — 



0.10 



0.24 
0.12 
0.14 



0.37 — 



2i/,, 



0.21 0.54 



0.23 
0.21 



c -I'np 
c +1'np 

-2!^NP 

+ 2j/„„ 



0.22 



0.27 
0.12 



0.24 0.26 0.21 0.12 



— — — — — 0.63 



1 1 

0.63 0.85 

0.40 0.68 

— 0.30 

— 0.17 



— — — 0.88 



V, 



K 



- Vr 



V^ + Vr 

V^ + Vr 

!^K + ^r - ^^NP 
V^+Vr + 2i/„p 



NP 



I ^NP 

-2i/,, 



0.89 0.86 
0.14 0.46 
— 0.15 



0.47 0.40 



0.39 
1 



0.26 0.45 

0.12 0.36 

— 0.25 



0.12 



0.27 



v^ + 2j^,. 

+ 2i/r - !^NP 
+ 2Ur + i^NP 
+ 2Ur - 2!/„„ 



1 1 

0.12 0.51 
— 0.19 



0.72 0.53 0.50 0.24 



0.22 



^The total power in each peak is normalised to the total power in the strongest peak produced by each effect. ■^Only 
peaks with total power > 10% of the strongest peak at frequencies < 2, 000 Hz are listed. 3The six columns for the HEG 
and MEG models list the power ratios for the prominent frequencies produced by the six different effects discussed in 
the text. These are (a) occultation of a uniformly bright neutron star surface; (b) occultation of a belt with a Gaussian 
brightness distribution of semi-width O.liJc around the stellar equator; (c) reflection of radiation from a uniformly bright 
neutron star surface; (d) radiation by a clump with a radius-independent proper luminosity; (e) radiation by a clump 
with a radius-dependent proper luminosity; (f) radiation by a clump in a tilted orbit that emits only while it is passing 
through a rarified disc of height h(T) = 0.05 r in the equatorial plane of the system. 

The actual power spectrum produced by a particular clump model is a linear combination of two or more of these 
spectra. For example, a clump that reflects radiation from the stellar surface [effect (e)] or produces a flare of emission 
as it passes through the orbital plane [effect (f )] will generally also occult the star [effects (a) and (b)] . 



at !/^p and i/^p — 'iv^^-p , but they have ~4-5 times less power 
than the peak at u^ + v,. + 2!/j^p and are no stronger than 
the peaks at Uj^ + Vr and v^ + Ur ~ 2//jjp. There is also a 
peak at v^ , but it has ten times less power than the peak at 
^K + ^r + 2!/pjp . The peak at 21/^^ has 12 times less power 
than the peak at u^+Ur + 2i'„p and 6 times less power than 
the peak at u^^ +2i'„p (it is too weak to be listed in Table M). 
There is no significant power at v^p. 

We emphasise that the actual power spectrum produced 
by a particular clump model is a linear combination of two or 
more of the spectra shown in Fig. [Ta For example, a clump 
that reflects radiation from the stellar surface [effect (c)] or 
produces a flare of emission as it passes through the orbital 
plane [effect (f)] will generally also occult the star [effects (a) 
and (b)]; indeed, occultation is typically the dominant effect. 



4.2 Nearly circular geodesies 

For completeness, we also computed the waveforms that 
would be produced by each of the six effects just discussed 
if the clumps were moving on nearly circular geodesies, even 
though such geodesies are excluded by the observed QPO 
frequency relations. Fig. hq shows power spectra of the typ- 
ical waveforms produced by such geodesies. The frequencies 
and relative total powers of the most prominent peaks are 
again listed in Table W. 

For orbits with eccentricities e = (ra — rp)/{rii + 
rp) < 0.1, there is (almost) always a strong peak at i/^, . Usu- 
ally, there is also a significant peak at i/j^ + Vi- Occultation 
can produce relatively strong peaks at i/„p and 2i/„p and 
significant peaks at the nodal sideband frequencies of v^ 
and i/j^ +1^1, but most effects do not produce peaks at these 
frequencies. 

We emphasise that for infinitesimally or only moder- 
ately eccentric orbits, r must be < {Re + d)/ cosi for the 



© 2000 RAS, MNRAS 000, 000-000 



Eccentric orbits and QPOs in neutron star X-ray binaries 21 



star to be even partially occulted. Hence, in order to pro- 
duce any occultation, such geodesies must have azimuthal 
frequencies 



10 km \ 3/2 



., >1550Hz(%/2cosi)3/2fi^') 



M 
2 Ma 



1/2 



(3) 



This requirement is inconsistent with the interpretation of 
the upper kilohertz QPO frequency as an azimuthal fre- 
quency, unless all the kilohertz QPO systems are highly 
inclined, which is unlikely and contradicts other evidence 
concerning their inclinations. 



5 GAS DYNAMICAL CONSTRAINTS 

The geodesic precession model assumes that at a particu- 
lar radius the accretion flow onto the neutron star changes 
from a disc flow to gas clumps moving on geodesies with the 
same or very similar apastron and periastron radii r^, and 
periastron Vp. It assumes further that the geodesies that are 
populated by clumps at any given time have the azimuthal, 
apsidal precession, and nodal precession frequencies required 
to explain the frequencies of the QPOs observed in the atoll 
and Z sources. Finally, the hypothesis assumes that as the 
accretion rate varies, the orbital parameters of the geodesies 
that are populated change in the way needed to produce the 
observed correlated variation of v^ or !^hbo, i^i, and V2- 

As we now discuss, the behaviour of the accreting gas 
required by the geodesic precession hypothesis is very dif- 
ficult to achieve in any physically consistent gas dynamical 
picture. For the purposes of the present section we set aside 
the results of the previous section, which showed that the 
frequencies generated by orbiting clumps are not the fre- 
quencies assumed by previous workers, and focus on the gas 
dynamical requirements of the model. 



5.1 Circularity of the flow in the disc 

We begin by showing that the stream lines in the accretion 
disc are necessarily nearly circular, in contrast to the highly 
eccentric geodesies required to produce orbital frequencies 
qualitatively similar to the observed frequencies of the kilo- 
hertz QPOs. The reason is that the apsidal precession rate 
varies rapidly with the radius of the geodesic. Hence any ec- 
centric, adjacent streamlines would intersect in less than an 
orbital period, creating strong shear stresses, pressure gra- 
dients, and possibly shock waves that will quickly circularise 
the motion. 

To see this, consider gas streaming along two nearby, ec- 
centric geodesies with periastron radii rpi and rp2 = Vpi + 
Svp. For the present purpose we assume the geodesies are 
ellipses in Euclidean space and describe each of them using 
the familiar radius-angle relation r = rp(l-|-e)/(l-|-ecos(^), 
where the angle is measured from the geodesic's peri- 
astron. Suppose the geodesies are initially concentric and 
aligned and have the same, finite eccentricity e. In the ab- 
sence of shear or pressure stresses, the different precession 
rates of the two geodesies will bring them into contact at 
the 'osculation' point after a time St. The angle between the 
axes of the two ellipses at this time can be obtained from the 
osculation point conditions ri — r2 and [dr/d(j>]i — [dr/d(j>]2, 



where the subscripts indicate the radius-angle relations for 
the two ellipses. 

Solving for 54> gives, to the lowest order in Srp/rp, 
smS(f> — {5rp/rp)(l — e'^Y'^ /e. The differential rate of peri- 
astron advance is Si/j^^ ~ {dv^p/drp)5rp. Using the lowest- 
order post-Newtonian expression for i^^p (see Markovic 
2000) we obtain 



5t = 



{54>/2^) 



1 (l+e)3/"(l-e) 



U/2 



[dUf^p/drp)5rp ISvr 



i^K,N 



(4) 



where i/K,N = (1 - ef^{GKlJrlf^^ /2tx is the Newtonian 
Kepler frequency. Equation M) shows that the two stream- 
lines will intersect in less than one orbital period even if 
they are only very slightly eccentric (e>0.02). Moderately 
or highly eccentric streamlines wil intersect in much less 
than one orbital period. 

In reality, the two streamlines will not cross. Instead, 
shear stresses and gas pressure forces will grow to oppose 
their differential apsidal precession, forcing the two stream- 
lines to deviate from purely geodesic motion. The result will 
be to circularise both streamlines within a time ~ 8t, i.e., in 
less than one orbital period. 



5.2 Effects of tidal forces 

The high coherence (Q ~ 100) of the kilohertz QPOs re- 
quires that the postulated gas clumps survive for at least 
Q orbital periods. However, the strong gravitational tidal 
force near the star will tidally disrupt clumps unless they 
are either confined by external forces or are very small. 

The rapid damping of noncircular flow in the inner disc 
makes persistence of the highly eccentric geodesic motion 
required by the geodesic precession model impossible un- 
less the clumps are isolated from the general gas flow, i.e., 
unless there is an extensive region around the star where 
the density of the interclump gas is negligible. The clumps 
therefore cannot be conflned by the surrounding gas. Nor 
can they be conflned by the stellar magnetic field, because 
the torque produced by a field strong enough to confine the 
gas in the clumps would be strong enough to cause them to 
deviate from purely geodesic motion. Therefore, the clumps 
can persist for Q orbital periods only if they are very small. 

To estimate how long a clump of radius d can per- 
sist, consider two elements of the clump moving on concen- 
tric orbits with the same eccentricity and orientation but 
slightly different periastra rp and r-p -\- 2d. Neglecting ap- 
sidal precession and other relativistic effects, the difference 
between the azimuthal frequencies of the two elements is 
5vk,k ~ (2d/rp)(3/2)i/K,N. Hence, after a time t the clump 
will be sheared out a distance Ad « rp2-K5vK,^t. The X- 
ray modulation produced by such a clump will certainly be 
markedly changed by the time the clump has been spread a 
quarter of the way around the star (i.e., over a phase inter- 
val of 7r/2). This will happen in less than Q orbital periods 
unless 



d 



< 



120 



8x 10" 



100 



(5) 



In the last expression on the right we have scaled Q by 
100, the minimum value required for consistency with the 
observed coherence of the kilohertz QPOs. Inequality (0) is a 
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Figure 17. Schematic depiction of the breakup of the accre- 
tion disc into dense clumps on nearly identical highly eccentric 
geodesies as required by the geodesic precession model. Here r^. 
and Tp are the apastron and periastron radii of the geodesic, risco 
is the radius of the innermost stable circular orbit, and the dark 
disc shows the size of a typical neutron star. In reality, geodesies 
near the star are very open rosettes rather than ellipses: a clump 
on these geodesies would go around the star 3—4 times as it moves 
from Ta, to rp . The figure shows that the clumps have relative ve- 
locities comparable to their azimuthal velocities, even if they are 
on geodesies with similar values of r^ and rp. 



strong upper bound on the size of a roughly spherical clump 
that would survive tidal disruption for Q orbital periods. As 
we show below, this bound is a very significant constraint 
on precessing clump models. 



5.3 Constraints on the gas density in the clumps 

The requirement that the interaction of clumps with the in- 
terclump gas not affect significantly their orbits for Q orbital 
periods leads to a lower bound on the density pc of the gas 
in the clumps relative to the density pi of the interclump 
medium. As shown above, the interclump gas must follow 
nearly circular orbits; hence any clumps on highly eccen- 
tric geodesies must move at supersonic speeds through the 
interclump gas. Assuming that it can be held together by 
some means, a gas clump of size d moving with supersonic 
velocity Vc — Vg through interclump gas of density pi will 
experience a drag force 



7^D=CDrf'pi(Vc-Vi) 



' Cod piV^ 



(6) 



where Co (~1) is the ballistic coefficient. A clump of mean 
density pc has momentum d^pcVc and hence persistence of 
geodesic motion for Q orbital periods would require a density 
ratio 



Pc rQCp 
Pi d 



> 10'' Co 



Viooy 



where in writing the last inequality on the right we have 
used condition (H). 

The requirement that the clumps not collide too fre- 
quently leads to an even more stringent lower bound on 
the density of the gas in the clumps relative to the mean 
density near the star. The mean clump-clump collision rate 
near the star is R = ricO-Vrei ~ ric^nd'^Vrei, where ric is the 
density of clumps near the star and v^ei is their typical rela- 
tive velocity. Even if the clumps all move on geodesies with 
the same values of rp and r-a, their periastron phases will 
be different (see Fig. y_7|) and hence they will approach each 
other with relative velocities that are comparable to vk ■ Now 
Uc = 3pc/{4:Tvpcd^), where pc is the mean mass density in the 
clumped region. Hence R ~ 3{pc/pc){vK/d). The observed 
coherence of the kilohertz QPOs requires 2-KQRrp/vK ^ 1, 
or 



pc 



^^»-'»"(i^r 



(8) 



where in the last inequality we have again used condition (m . 
These dynamical constraints show that the density con- 
trast in the clumps required by the geodesic precession hy- 
pothesis is enormously greater than the density contrasts 
and magnetic field fluctuations considered in the magneto- 
spheric (Alpar & Shaham 1985; Lamb et al. 1985) and sonic- 
point (Miller et al. 1998) beat-frequency models of QPOs; 
in the latter models, the fluctuations have relative ampli- 
tudes ^ 1 and move with the mean flow. How clumps with 
such extremely high densities would form or persist in a gas 
dynamical flow and how they would be replenished are un- 
known. A large number of clumps would have to be produced 
to carry a substantial fraction of the accretion flow and pro- 
duce significant X-ray modulation. Clumps would have to be 
formed from the gas in the region where the orbits are nearly 
circular and then "shot" quickly into the required highly ec- 
centric orbits in an inner, highly rarefied region near the star, 
in order to avoid interaction with the gas moving on nearly 
circular orbits. However, any forces capable of producing the 
required sudden deviation from circular motion would also 
be strong enough to disrupt the clumps and cause them to 
deviate strongly from purely geodesic motion. As noted ear- 
lier, what picks out the particular geodesies required by the 
geodesic precession hypothesis is also unknown. 



5.4 Implications for QPO amplitudes 

The small size necessary for clumps to persist as long as 
required by the observed coherence of the kilohertz QPOs 
means that the projected area of a clump is only a very 
small fraction of the projected area of the stellar surface. 
The fractional area of a clump is 



/ 



TT<f_ 



< 



d^ 



< 7 X 10" 



(9) 



This is also an upper bound on the fractional modulation 
of the X-ray fiux that can be produced by clumps, whether 
the modulation is caused by occultation of the surface of the 
star, refiection of radiation from the star, or emission from a 
clump at about the same temperature as the stellar surface. 



(7) 
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6 DISCUSSION AND CONCLUSIONS 



If it could be established that the frequencies of the most 
prominent humps, bumps, and QPOs observed in power 
spectra of the X-ray emission from neutron star and black 
hole LMXBs are simply the characteristic frequencies of gas 
clumps orbiting these objects and that hydrodynamic, ra- 
diation, and magnetic forces have no significant effects, so 
that their motion is purely geodesic motion, this would be 
an extraordinarily important result. In particular, it would 
mean that measurements of these frequencies would provide 
a very simple way to explore the effects of strong-field grav- 
ity and the properties of compact objects. This is why we 
have carefully investigated this suggestion. Here we discuss 
our principal results and summarize our conclusions. 



6.1 Geodesies and frequency relations 

We first explored whether sequences of geodesies could be 
found with azimuthal and apsidal precession frequencies z/^p 
and i/j^ that agree, respectively, with the observed frequen- 
cies vi and V2 of the lower and upper kilohertz QPOs. We 
also investigated whether geodesic sequences could be found 
for which the first overtone 21^^^^ of the nodal precession 
frequency and u^ and Uj^p agree, simultaneously, with the 
frequencies of the low-frequency and kilohertz QPOs. Here 
we briefly summarise the results of these investigations. 



6.1.1 Infinitesimally and moderately eccentric geodesies 



In Section i.l we explored whether it is possible to con- 
struct sequences of infinitesimally or moderately eccentric 



geodesies that give u^ 



relations consistent with the ob- 



served correlations between ui and 1/2- 

We first investigated infinitesimally eccentric geodesies. 
The geodesic precession hypothesis makes no predictions 
about the apastron and periastron radii ra and rp of the 
geodesies or how they vary, except that they must (obvi- 
ously) be larger than the equatorial radius of the star. In 
exploring infinitesimally eccentric geodesies, we set rp = ra 
and treated ra as a free parameter, subject only to the re- 
quirement that it exceed the radius of the innermost stable 
circular orbit riaco. The parameters in this model are there- 
fore the EOS, mass, and spin frequency of the star, and 
an orbital radius for each pair of kilohertz QPO frequencies 
that was fit. 

We found that the ^^ap'^^k relations given by the best- 
fitting sequences of infinitesimally eccentric geodesies are 
much steeper than the observed Au = 1/2 — 1^1 vs. 1^2 corre- 
lation. Our i^AP'^^K relations are consistent with the Uj^p-v^^ 
relations reported by SV, which are also much steeper than 
the frequency data presented in their paper. The qualitative 
disagreement between the best-fitting frequency relations for 
infintesimally eccentric geodesies and the kilohertz QPO fre- 
quency data is indicated by the very large values of x^/dof 
found ( ^ 35 for the Z source Sco X-1 and > 10 for the atoll 
source 4U 1728-34). 

We note that the qualitative disagreement between the 
frequency relations for infinitesimally eccentric geodesies 
and the frequency data is a serious difficulty for the res- 
onant ring model recently proposed by Psaltis & Norman 
(2000) to explain the QPOs in both neutron star and black 



hole LMXBs. Not only does this disc model predict a dense 
spectrum of power spectral peaks, in contrast to the small 
number of peaks observed, the frequencies are necessarily 
the frequencies of nearly circular geodesies. In order for 
"hydrodynamic corrections" (i.e., gas pressure forces) to be 
large enough to shift the frequencies by ~15% to bring them 
into agreement with the data, the disc height would have to 
be ~ 0.2 r, invalidating the narrow ring approximation on 
which the model is based. The behaviour of the disc fre- 
quencies would then differ qualitatively from the behaviour 
predicted by the narrow ring model. 

Next, we investigated geodesies with finite but moder- 
ate eccentricities. Once geodesies with finite eccentricities 
are considered, the geodesic precession hypothesis does not 
by itself predict a relation between i/^p and v^ (see the dis- 
cussion in Section nl). In order to create a relation between 
Vj^p and i/j^ , SV proposed that some physical mechanism 
keeps rp constant in a given source while ra varies. They 
suggested that rp may remain constant because it is equal 
to the radius riaco of the innermost stable circular orbit, the 
radius r™ at which the magnetic field of the star first couples 
to the accreting gas, or the equatorial radius R^ of the neu- 
tron star. However, there is no physical reason why the peri- 
astron radius of an eccentric geodesic cannot be significantly 
smaller than risco (see Markovic 2000). The magnetic cou- 
pling radius decreases as the mass accretion rate increases 
(Lamb 1989; Ghosh & Lamb 1992; Miller et al. 1998; Psahis 
et al. 1999a). Moreover, if the stellar magnetic field is strong 
enough to truncate orbits at r^, the motion of the accreting 
gas will not be purely geodesic motion for a substantial dis- 
tance outside rm. The equatorial radius Re provides a lower 
bound on rp, but otherwise provides no constraint. 

In exploring sequences of geodesies with finite eccen- 
tricities, we simply followed SV, who assumed that rp is 
constant in a given source but that it can be chosen freely 
for each source to give the best possible fit of the i^ap'^k 
and 2i^^p-!/j^ relations to the observed 1/1-^2 and i^^^-i^^ cor- 
relations. We also assumed that for each source ra can be 
varied freely to give the best possible agreement between 
the allowed frequency relation and the observed frequency 
correlation, requiring only that it exceed rp. Therefore the 
parameters in this model are the EOS, mass, and spin fre- 
quency of the neutron star, the periastron radius rp, and an 
apastron radius for each pair of kilohertz QPO frequencies 
that was fit. 

We found that the Vj^p-v^^ relations given by the best- 
fitting sequences of moderately eccentric geodesies, defined 
as geodesies with r^jr^ < 1.5, also disagree qualitatively 
with the kilohertz QPO frequency data. As we explained 
in detail in Section BJ| and the Appendix, SV came to a 
different conclusion in their paper proposing the geodesic 
precession model because they computed the azimuthal fre- 
quency v^ incorrectly. When computed correctly, the az- 
imuthal frequencies of the geodesies cited by SV as giving 
agreement with the Sco X-1 frequency data miss most of the 
data points by ~100Hz. Our investigation shows that the 
best-fitting sequences of moderately eccentric geodesies dis- 
agree qualitatively not only with the kilohertz QPO frequency 
data on Sco X-1 but also with the data on other sources, re- 
gardless of the neutron star model and spin rate assumed. 
The reason sequences of such geodesies are unable to fit the 
data is that although the frequency difference Av = 1/2 — 1^1 



© 2000 RAS, MNRAS 000, 000-000 



24 D. Markovic and F. K. Lamb 



decreases with increasing 1/2, the radial epicyclic frequency 
Ur = Vj^ — Vj^p given by such geodesies decreases much more 
rapidly with increasing u^ for any realistic neutron star and 
any allowed spin rate. 



6.1.2 Highly eccentric geodesies 



In Section 3.2 we explored whether dropping the require- 
ment that geodesies be only moderately eccentric would al- 
low acceptable fits to the kilohertz QPO frequency data. We 
were able to find sequences of highly eccentric geodesies that 
give i^Ap'^K relations qualitatively similar to the frequency 
data on some sources, but only by using geodesies with very 
high eccentricities {ra/rp > 3). We were able to find geodesic 
sequences that agree qualitatively with the frequencies ob- 
served in Sco X-1 and 4U 1728-34, provided that the EOS 
of neutron star matter is similar to the UU or A18+\JlX+5vb 
equations of state, although the fits are not formally ac- 
ceptable. In general, frequency relations similar to the ob- 
served frequency correlations are possible only if values of 
Tp smaller than risco are allowed. We were able to find se- 
quences of highly eccentric geodesies that give formally ac- 
ceptable fits to the GX 17+2 and GX 340-f kilohertz QPO 
frequencies; however, only a few frequency measurements 
have been reported for these sources and the uncertainties 
are large. The best-fitting geodesic sequence appears to be 



inconsistent with the GX 5—1 frequency data. In Section 3.3 
we showed that geodesic frequencies other than j/^p and i/^ 
are inconsistent with the frequencies of the kilohertz QPOs. 
In order to test the suggestion of SV and SVM that the 
low-frequency QPOs in the atoll sources and the HBOs in 
the Z sources are the first overtones of nodal precession fre- 
quencies, in Section 3.4 we investigated whether the nodal 
precession frequencies of the sequences of best-fitting highly 
eccentric geodesies found in Section Bj are consistent with 
the frequency behaviour of the low-frequency QPOs. We 



found that the predicted u^ 



and 2u,, 



relations are 



qualitatively different from the i'l-i'k ^^'^ vh^q-v^ correla- 
tions observed in 4U 1728-34, GX 5-1, and Sco X-1. The 



predicted 2v^ 



relation is somewhat similar to the fre- 



quency correlation observed in GX 17-1-2, if the spin rate of 
this neutron star is 660 Hz. However, the value of x^/dof 
is large (4.3) and the required stellar mass is equal to the 
maximum stable mass for this spin rate. 

we considered the unusual source 



Finally, in Section 3.5 



Cir X-1, which is thought to be a neutron star (see van der 
Klis 2000). The frequencies v-^ and v-^^ of its low- frequency 
QPO and high frequency 'hump' both vary by more than 
an order of magnitude (their reported ranges are 1-12 Hz 
and 20-200 Hz, respectively; see Tennant, Fabian & Shafer 
1986; Tennant 1987; Shirey 1998; Shirey et al. 1996, 1998, 
2000). This makes Cir X-1 a crucial testing ground for any 
model that seeks to explain the neutron star and black hole 
QPOs within a single framework. Indeed, without the Cir X- 
1 data, the evidence for a connection between the black hole 
and neutron star QPOs would be much weaker (see Fig. 2 
of Psahis et al. 1999b). 

We first investigated whether a sequence of moderately 
eccentric geodesies can be constructed that gives a 2u-^-p-Vj^-p 
relation similar to the Vi^-v^ correlation observed in Cir X-1 
(I'p^p is much too low to be consistent with v-^). We find that 
sequences of moderately eccentric geodesies around rapidly 



spinning {u^ > 650 Hz) neutron stars can give 2Vj^p-u^-p re- 
lations that are roughly consistent with the observed u-^- 
i/jj correlation for v-^ < 100 Hz. However, only moderately 
eccentric geodesies around very low mass {M <,1.3 Mq for 
EOS UU) neutron stars spinning near breakup give rela- 
tions roughly consistent with both the lower-frequency data 
and the //^ ~ 12 Hz, v^^ « 200 Hz measurement by Tennant 
(1987). No quantitative comparison is possible, because no 
errors have been reported for most of the frequency mea- 
surements. We then explored whether using sequences of 
highly eccentric geodesies would make it easier to construct 
a frequency relation similar to that observed in Cir X-1. 
We found that the situation for highly eccentric geodesies 
is similar to that for moderately eccentric geodesies. The 
lower-frequency data favour moderately eccentric geodesies 
that are extremely eccentric (ra/rp ^ 6) and far from the 
neutron star (ra/-Rc ^ 10). Inclusion of the higher-frequency 
measurement of Tenant (1987) again makes it much more 
difficult to find a sequence of geodesies that gives a frequency 
relation similar to the observed correlation. Only highly ec- 
centric geodesies around stars with very low masses and very 
high spin rates give frequency relations that pass near all the 
measurements. 

In the majority of cases, the values of rp required to pro- 
duce Vj^p-i'y^ relations similar to the observed kilohertz QPO 
frequency corelations are substantially larger than the equa- 
torial radius Re of the corresponding neutron star model (for 
example, in GX 340-1-0 the required value of rp is more than 
2Re). Where a marginally bound geodesic exists, rp is also 
significantly larger than rmb for the best-fit HEG sequences. 
Hence the geodesic precession hypothesis does not provide 
any obvious explanation of why the periastron radius should 
have the required value or of why it should remain fixed in a 
given source. As discussed in Section y, it seems impossible 
to overcome this difficulty for any realistic EOS. 

The situation for Cir X-1 is even more serious, because 
geodesies with apastron radii ra ~ 5-50 Re are required to fit 
the frequency data. Clumps so far from the star are unlikely 
to produce significant X-ray modulation. 



6.2 Expected frequencies and power spectra 

Up to this point in our analysis we simply assumed that 
the dominant modulation frequencies produced by clumps 
moving on geodesies around a neutron star would be i/^, 
J^j^p, and 2^- p, as proposed by SV, SVM, and Karas (1999). 
In Section ]d we carried out extensive numerical simulations 
to determine the waveforms and power spectra that would 
actually be produced by geodesic motion of clumps. We con- 
sidered six effects: occultation of the whole stellar surface, 
occultation of a bright equatorial band, reflection of radi- 
ation by clumps, two models of emission of radiation by 
clumps, and interaction of clumps with a rarified gaseous 
disc in the rotation equator of the star. Any orbiting clump 
model will generally produce a power spectrum that is a su- 
perposition of the power spectra produced by these effects. 
We find that, contrary to what was simply assumed 
previously, f^^, Vj^p, and 2z/pjp are not the most prominent 
frequencies generated by orbiting clumps. Instead, most of 
the power generated by such clumps is typically in peaks of 
roughly equal power at the radial epicyclic frequency Vr, at 
i^K + '^r, and at u^ + 2ur (see Table 41). As explained in Sec- 
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tion kl the reason for this is that for the general relativistic 
geodesies that are relevant, the only strictly periodic motion 
is the radial epicyclic motion. In particular, the azimuthal 
motion is not periodic. There is usually power at u^, but 
it is always < 25% of the power in each of the dominant 
peaks and there are often many other peaks as strong or 
stronger. For relevant geodesies, the power at !/^p and at 
I'pjp produced by most mechanisms is < 10% of the power 
in the dominant peaks. We were unable to find any effect 
or geometry that produces a significant peak at 2iy^p . Thus, 
the frequencies of the significant peaks in the power spectra 
that would be generated by orbiting clumps are qualitatively 
different from the frequencies of the peaks observed in power 
spectra of the kilohertz QPO sources. 

Periodic modulation of the X-ray emission by orbiting 
clumps is possible only if the inclination i of the system, 
which is expected to be the same as the inclination of the 
star's spin axis and the normal to the accretion disc, is suf- 
ficiently large. More generally, the frequency structure and 
the relative strengths of the oscillations at different frequen- 
cies produced by this mechanism depend sensitively on the 
inclination of the system. However, there is no evidence for 
any inclination-dependence of the frequency pattern in the 
low- or high-frequency QPOs, or even in their relative am- 
plitudes. This is difficult to understand in orbiting clump 
models. 



6.3 Dynamical considerations 

In Section o we investigated the dynamical constraints on 
orbiting clump models and the X-ray modulation generated 
by such models. We find that in order for clumps to persist 
in geodesic motion for ~100 orbits, as required by the coher- 
ence of the kilohertz QPOs, their radii must be <^ 10~^ of 
the periastron radius. In order to avoid rapid inspiral caused 
by drag forces, the density of the gas in the clumps must be 
^ 10^ times the density of the interclump medium. In order 
for clumps to avoid colliding with each other too frequently, 
the density of the gas in the clumps must be S> 10^ times 
the mean density near the star. These densities are extreme. 
As discussed in Section BI the small size required for 
clumps to persist as long as required by the observed coher- 
ence of the kilohertz QPOs means that the fractional mod- 
ulation of the X-ray flux that can be produced by clumps 
is <IC 10~^. This upper bound is many orders of magnitude 
smaller than the observed fractional amplitudes, which are 



-10" 



-10~^ (vander Klis 2000). 



The very low interclump gas density and the small cov- 
ering factor by clumps that is required by the geodesic pre- 
cession hypothesis means that radiation can escape from the 
clumps and from the stellar surface without suffering signifi- 
cant scattering or absorption. Hence scattering or absorption 
can have little effect on the frequency content of the wave- 
forms produced by the clumps. Therefore all the peaks at 
all the frequencies we have discussed should be observable. 
However, most have not been detected (see, e.g., Jonker et 
al. 1998). 

A related point is that even weakly nonuniform emission 
from the stellar surface should produce readily observable 
periodic X-ray oscillations at the stellar spin rate. Again, 
such oscillations have not been detected (see van der Klis 
2000). 
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A major deficiency of the geodesic precession model is 
that no attempt has been made to explain what causes the 
gas in the accretion disc, which is moving along nearly cir- 
cular streamlines, to break up into clumps moving on highly 
eccentric geodesies. This would require the sudden removal 
of a substantial fraction of the angular momentum of the 
gas in the clumps. Nor has there been any attempt to ex- 
plain why and how clumps are formed only on geodesies 
with nearly the same values of r^ and Vp, why these spe- 
cial geodesies have similar frequencies and obey similar fre- 
quency relations in sources with very different properties 
(e.g., accretion rates that differ by factors ~30), or why the 
properties of the special geodesies change with the accretion 
rate in just the way needed to explain the observed changes 
in the frequencies of the QPOs. Finally, no explanation has 
been given of how the clumps dissipate and how the accret- 
ing gas reaches the stellar surface (injection of clumps onto 
geodesies that do not intersect the stellar surface does not 
by itself produce any accretion onto the star) . 



6.4 Conclusions 

We conclude that there are significant difficulties in inter- 
preting the kilohertz and low-frequency QPOs in the atoll 
and Z sources as consequences of gas clumps moving along 
geodesies. 

The best-fitting geodesic sequences that can be con- 
structed using moderately eccentric geodesies give frequency 
relations that differ qualitatively from the kilohertz QPO 
frequency data. Frequency relations qualitatively similar 
to the kilohertz QPO data on some sources can be con- 
structed, but only if highly eccentric geodesies are used. Such 
geodesies are very difficult to accommodate in any physi- 
cally consistent gas dynamical picture of the accretion flow 
near the star. The best-fitting 2u^p-iy^ and Vj^p-v^^ relations 
disagree qualitatively with the low-frequency and kilohertz 
QPO frequency data, for any realistic neutron star models 
and allowed spin frequencies. Other geodesic frequencies are 
qualitatively inconsistent with the observed frequencies. 

Modulation by orbiting clumps would not produce 
strong peaks in the power spectra at the frequencies i^^p , v^ , 
and 2!/j^p as required by the geodesic precession hypothesis. 
Instead, the strongest peaks would be at the radial epicyclic 
frequency v^ and at u^^ + Vi and i/j^ -I- 2vi , frequencies that 
are inconsistent with the frequencies observed. 

Finally, the dynamical constraints on orbiting clump 
models make it very difficult to see how clumps of the re- 
quired density could be formed on highly eccentric geodesies, 
or how X-ray modulation with the amplitudes observed 
could be generated by orbiting clumps. 
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APPENDIX A: KEPLERIAN (AZIMUTHAL) 
FREQUENCIES OF ECCENTRIC GEODESICS 

In order to clarify the difference between the definition of the 
Keplerian (azimuthal) frequency of eccentric geodesies and 
the quantity used by Stella & Vietri (2000; hereafter SV), 
we first discuss the definition of the Keplerian frequency, 
then derive expressions for the radial epicyclic and Keple- 
rian frequencies of eccentric geodesies in the Schwarzschild 
spacetime, and finally explain how we were able to repro- 
duce the curve shown by SV in their Fig. 2. The expressions 
needed are given by Markovic (2000). However, in order to 
facilitate comparison with the expressions given by SV, in 
this appendix we use the notation of SV, whose derivation 
follows closely the discussion in Chandrasekhar (1983, pp. 
100-105). 



Al Definition of the Keplerian frequency 

For a test particle moving on an eccentric geodesic, the rate 
of azimuthal phase advance varies with its radius, and hence 
only the average rate of azimuthal phase advance is mean- 
ingful. This is the quantity conventionally denoted 27ri^^ or 
27r!/j^. In the Schwarzschild spacetime, eccentric geodesies 
are not closed (because of apsidal precession), the azimuthal 
motion of a particle moving on such a goedesic is not peri- 
odic, and the average rate of azimuthal phase advance is not 
proportional to the inverse of the time required for the par- 
ticle to traverse a 27r interval in azimuthal phase (see, e.g.. 
Section 14 of Landau & Lifshitz 1976) . However, the radial 
motion is periodic, with period l/t'r, where Ur is the radial 
epicyclic frequency. The rate at which a particle's azimuthal 
phase advances depends only on its radius and hence is also 
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periodic with period l/vr- The average rate at which a par- 
ticle's azimuthal phase advances can therefore be computed 
by calculating the phase advance during any time interval 
that is a multiple of half the radial epicyclic period and then 
dividing by the time interval. For example, if /S.cj> is the az- 
imuthal phase advance during each radial epicyclic period, 
then Vijy = {A(f)/2n) Vr- 

Instead of using the usual definition of v^,, SV appear 
to have assumed that i/^ is the inverse of the time required 
for the azimuthal phase to advance by 27r. However, as just 
explained, the rate of azimuthal phase advance of a particle 
moving on an eccentric geodesic in the Schwarzschild space- 
time depends on the instantaneous radius of the particle 
and the azimuthal motion is not periodic. Hence the time 
required for the particle's azimuthal phase to increase by 2tt 
depends on the assumed initial position of the particle and 
is not even well-defined; it is not the inverse of the Keplerian 
frequency. As Fig. Al shows, the time required for the az- 
imuthal phase to advance by 2n depends sensitively on the 
starting position used in computing it, even when the eccen- 
tricity is relatively small. The non-uniqueness of this time 
is a consequence of the nonperiodic nature of the azimuthal 
motion. (The nonperiodic nature of the azimuthal motion 
is also the reason why the peak at u^ in power spectra of 
the X-ray flux modulation produced by eccentric geodesic 
motion is not very prominent.) 



A2 Calculation of the Keplerian and radial 
epicyclic frequencies i/^ and i/r 

For a non-rotating {j — 0) spherical star, equations (20)- 
(22) of Markovic (2000), which describe the geodesic mo- 
tion of a particle of unit proper mass, simplify considerably. 
Following Chandrasekhar (1983), we use the radial variable 
u = 1/r, where r is the radius in Schwarzschild coordi- 
nates. Then (see Chandrasekhar 1983, Section 19; Shapiro 
& Teukolsky 1983, Ch. 12) 

ii'^ = (2L^u^ - L'^u^ +2u + E'^ -l) m" 
E 
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where I = 2/(ui + U2) and e = (ra — rp)/{r^ + rp) = {u2 — 
ui)/{u2 -I- Ml) is the ecce nt ricit y. 

From equations ( |Al[ )-( p\^ , we obtain 



2u 



= Lu 



(Al) 



where E and L are, respectively, the conserved energy and 
angular momentum of the particle. Here all distances are in 
units of GM/c^ and all t imes, in units of GMJC' . 
From equations (|A1| ) it follows that 



2(u — u\){u — U2)(U — U3), 



(A2) 



where u\ = 1/ra, U2 — i/rp, and us > U2 > ui. Introducing 
/x = l/l and evaluating the coe fficients of the polynomials 
in u in equations (Al) and (A2), one finds 



Ul = 

and 

J_ 
I2 
E^ 

X2 



U2 = 



1 + e 



1 2 

"^ ^ 2 - T 



(A3) 






-m(3 + 0] 

- 2^^f - 4e^f^^] 



(A4) 



— 1 = l-2/i(3 + ecosx) 

= (l-6M + 2eAt) [1-k2cos2(x/2)] , (A5) 

where k^ = 4e/i/(l — 6/i + 2e/i) and x is defined implicitly 
by 

_ 1 + e cos X 



I 



(A6) 



Equations (Al) and (AS) yield 

dt dt d(j) 

dx dcj> dx 



\l-2u) [lu^J 



[1 - 2^(3 + e cos x)]'/' 
Z3/2 [(1 - 2iif ^ Ae^y.'^] "'' (1 + ecosx)"' 
[1 - 2^(1 +£ cos x)][l - 2^(3 -^ ecosx)]i/2 



(A7) 



As X increases from to 2ii, the particle executes one 
full period of its radial motion, moving from its periastron 
radius rp = l/u2 to the apastron radius r^ = 1/ui and back 
to Tp again. Hence the radial epicyclic frequency Vi is given 
by 



1 



, .^1* 



(A8) 



One can relate Ur to the Newtonian frequency un ~ 
a~^'^/2Tv using the expression a = (l/wi -I- l/u2)/2 = 
1/(1 — e^) for the semi-major axis a of the geodesic. The 
result is 

a^'^N ,, \.,,. - _ ..^ ,, ./ . , (A9) 



where 
/(e,m) 



(l_e2)3/2[(i_2^)2-4eV]'^'/(^'/^) 

(l + ecosx)~^rfx 

[1- 2^(1 + ecosx)][l- 2^(3 + ecosx)]i/2 ' 

(AlO) 



As the particle executes one period of its radial epicyclic 
motion, its azimuth (p increases by (see eq. A5) 



A0: 



yi — 6/i + 2ey 



dx 



[1 



Hx/2)] 



1/2 



(All) 



The mean rate of advance of the particle's azimuth is 
27r!/0 = z^r Acjf). Hence 

l^^^Ur-^. (A12) 

Quite generally, A<j!> > 2tt, so the geodesies are not closed. 

A3 Calculation of SV 



Equations ( |A9|) and ( A.10 ) for the radial epicyclic frequency 
Vr are the same as equations (2) and (3) given by SV. How- 
ever, instead of com puti ng th e azimuthal frequency u^ given 
by equations (All) and ( A.12), SV appear to have computed 
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Figure Al. The data points in the left panel show the Av-i/2 correlation observed in ScoX-1. The solid curve labeled K in the left panel 
shows how the radial epicyclic frequency Vr varies with the Keplerian (azimuthal) frequency vy, for the sequence of geodesies specified 
in the caption of Fig. 2 of Stella & Vietri (1999). For these geodesies, the eccentricity e increases with increasing I'r, as shown in the 
right panel. The dotted curve labeled ■k in the left panel shows Vr as a function of l/t,r, the inverse of the time required for a particle to 
advance in azimuthal phase from tt radians before apastron to vr radians after apastron. The short dashed curve shows v,. as a function 
of l/isTTi the inverse of the time required to advance from vr radians after apastron to Svr radians after apastron. The long dashed curv 
shows Vr as a function of I/Jstt , the inverse of the time required to advance from Stt radians after apastron to Stt radians after apastron. 
The dash-dotted curve shows Vr as a function of l/t7,r , the inverse of the time required to advance from Stt radians after apastron to Ttt 
radians after apastron. 



the time required for a particle to advance through an az- 
imuthal interval of 2-k centered about the azimuthal phase 
of its apastron and then to have identified the inverse of this 
time with v^, which it is not. 

If a particle has X — Xo initially, its azimuthal phase will 
have increased by 27r at x = Xi > where xi is given implicitly 
by 



27r: 



\J\ — 6/i + 2e/i 



Xl 



dx 



Xo 



[1-k2cos2(x/2)] 



1/2 



(A13) 



In the Schwarzschild spacetime, the radial motion of a par- 
ticle is always slower than its azimuthal motion and hence 
Xi ~ Xo < 27r. As noted above, the rate at which the par- 
ticle's azimuthal phase advances depends on its radius and 
its motion in the azimuthal direction is nonperiodic. Hence, 
the time required for the azimuthal phase of a particle to 
advance by 27r depends on its initial position, i.e., on xo. 

According to SV, the azimuthal frequency is given by 
the inverse of their equation (2) (our eq. [|A9||) , but with 
/(e,/i) given by their equation (3) (our eq. |A1C| ]), the inte- 
gral in their equation (3) taken over to Xc rather than to 
27r, and Xc given implicitly by their eq. (4). This 'azimuthal 
frequency' is the inverse of the time required for the particle 
to traverse the azimuthal phase interval of 27r centered on 
the particle's apastron, if the integration in SV's eq. (3) is 
meant to be over vr — 2xc to tt -I- 2xc, rather than to Xc, 
and their expression for k is instead an expression for fc^.N 

T SV follow Chandrasekhar's (1983) notation, but use x instead 
of X- Hence they use the variable 7 = 7r/2 — x/2 rather than x 



Indeed, we were able to reproduce the v-^ vs. 'azimuthal fre- 
quency' curve shown in Fig. 2 of SV (see below) by using as 
the 'azimuthal frequency' the inverse of the time required 
for X to increase from Xo to Xi, where Xi ~ ^ = ■"" ~ Xo- (Xo 
is related to Xc by Xo = ■"" ~ 2xc.) This is the inverse of the 
time required for the particle to traverse the 2-k phase in- 
terval centered on its apastron phase (see eq. |A6]). Fig. Al 
shows that the times required for the azimuthal phase of the 
particle to advance by successive additional intervals of 2-k 
are all different from one another. 

The rate of azimuthal phase advance is smallest at apas- 
tron, and hence the azimuthal phase interval used by SV 
gives the smallest possible inverse time and thus the flattest 
possible i/r-inverse time curve. This flattest possible curve 
is shown as the dotted curve in Fig. Al, for the parameters 
Stella & Vietri list in the caption of their Fig. 2. Fig. Al 
shows that any other choice of inverse time would have given 
a i^r-inverse time curve much further from the Sco X-1 data. 
These inverse times are not relevant in any case, because 
none of them are equal to th e az imuthal frequency of the 
particle. As shown in Section 3.1, the best-fitting Vr-v^ re- 



lation is qualitatively different from the correlation between 
the kilohertz QPO frequencies observed in Sco X-1. 



in their eq. (4), which corresponds to our eq. (A13). Also, they 
write their eq. (4) in terms of the integral over only the first half 
of the 27r phase interval centered at the particle's apastron. 
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